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ANALYSIS OF LANDING-GEAR BEHAVIOR 1 


By Benjamin Milwitzky and Francis E. Cook 


SUMMARY 

This report presents a theoretical study of the behavior of the 
conventional type of oleo-pneumatic landing gear during the 
process of landing impact. The basic analysis is presented in 
a general form and treats the motions of the landing gear prior 
to and subsequent to the beginning of shock-strut deflection. In 
the analysis of the first phase of the impact the landing gear is 
treated as a single-degree-of -freedom system in order to deter- 
mine the conditions of motion at the instant of initial shock-strut 
deflection , after which instant the landing gear is considered as 
a system with two degrees of freedom. The equations for the 
two-degree-of-freedom system consider such factors as the 
hydraulic {velocity square) resistance of the orifice , the forces 
due to air compression and internal friction in the shock strut , 
the nonlinear force-deflection characteristics of the tire , the wing 
lift , the inclination of the landing gear , and the effects of wheel 
spin-up drag loads. 

The applicability of the analysis to actual landing gears has 
been investigated for the particular case of a vertical landing 
gear in the absence of drag loads by comparing calculated 
results with experimental drop-test data for impacts with and 
without tire bottoming. The calculated behavior of the landing 
gear was found to be in good agreement with the drop-test data. 

Studies have also been made to determine the effects of varia- 
tions in such parameters as the dynamic force-deflection 
characteristics of the tire , the orifice discharge coefficient, and the 
polytropic exponent for the air-compression process , which 
might not be known accurately in practical design problems. 

The study of the effects of variations in the tire characteristics 
indicates that in the case of a normal impact without tire 
bottoming reasonable variations in the force-deflection character- 
istics have only a relatively small effect on the calculated behavior 
of the landing gear. Approximating the rather complicated 
force-deflection characteristics of the actual tire by simplified 
exponential or linear-segment variations appears to be adequate 
for practical purposes. Tire hysteresis was found to be 
relatively unimportant. In the case of a severe impact involving 
tire bottoming, the use of simplified exponential and linear- 


segment approximations to the actual tire force-deflection 
characteristics, which neglect the effects of tire bottoming, 
although adequate up to the instant of bottoming, fails to indicate 
the pronounced increase in landing-gear load that results from 
bottoming of the tire. The use of exponential and linear- 
segment approximations to the tire characteristics which take 
into account the increased stiffness of the tire which results from 
bottoming, however , yields good results. 

The study of the importance of the discharge coefficient of the 
orifice indicates that the magnitude of the discharge coefficient 
has a marked effect on the calculated behavior of the landing 
gear; a, decrease in the discharge coefficient {or the product of the 
discharge coefficient and the net orifice area) results in an 
approximately proportional increase in the maximum upper- 
mass acceleration . 

The study of the importance of the air-compression process 
in the shock strut indicates that the air springing is of only 
minor significance throughout most of the impact and that 
variations in the effective polytropic exponent n between the 
isothermal value of 1 .0 and the near-adiabatic value of 1 .3 have 
only a secondary effect on the calculated behavior of the landing 
gear. Even the assumption of constant air pressure in the strut 
equal to the initial pressure, that is, n=0, yields fairly good 
results which may be adequate for many practical purposes. 

In addition to the more exact treatment, an investigation has 
been made to determine the extent to which the basic equations 
of motion can be simplified and still yield acceptable results. 
This study indicates that, for many practical purposes, the 
air-pressure force in the shock strut can be completely neglected, 
the tire force-deflection relationship can be assumed to be linear, 
and the lower or unsprung mass can be taken equal to zero. 
Generalization of the equations of motion for this simplified 
system shows that the behavior of the system is completely 
determined by the magnitude of one parameter, namely the 
dimensionless initial-velocity parameter. Solutions of these 
generalized, equations are presented in terms of dimensionless 
variables for a wide range of landing-gear and impact parameters 
which may be useful for rapidly estimating landing-gear 
performance in preliminary design. 


1 Supersedes NACA TN 2755, “Analysis of Landing-Gear Behavior” by Benjamin Milwitzky and Franeis E. Cook, 1952. 
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INTRODUCTION 

The shock-absorbing characteristics of airplane landing 
gears are normally developed largely by means of extensive 
trial-and-error drop testing. The desire to reduce the ex- 
pense and time required by such methods, as well as to pro- 
vide a more rational basis for the prediction of wheel-inertia 
drag loads and dynamic stresses in flexible airframes during 
landing, emphasizes the need for suitable theoretical methods 
for the analysis of landing-gear behavior. Such theoretical 
methods should find application in the design of landing 
gears and complete airplane structures by permitting 

(a) the determination of the behavior of a given landing- 
gear configuration under varying impact conditions (velocity 
at contact, weight, wing lift, etc.) 

(b) the development of a landing-gear configuration to 
obtain a specified behavior under given impact conditions 

(c) a more rational approach to the determination of wheel 
spin-up and spring-back loads which takes into account the 
shock-absorbing characteristics of the particular landing gear 
under consideration 

(d) improved determination of dynamic loads in flexible 
airplane structures during landing. This problem may be 
treated either by calculating the response of the elastic sys- 
tem to landing-gear forcing functions determined under the 
assumption that the airplane is a rigid body or by the simul- 
taneous solution of the equations of motion for the landing 
gear coupled with the equations representing the additional 
degrees of freedom of the structure. In many cases the former 
approach should be sufficiently accurate, but in some 
instances, particularly when the landing-gear attachment 
points experience large displacements relative to the nodal 
points of the flexible system, the latter approach, which takes 
into account the interaction between the deformations of the 
structure and the landing gear, may be required in order to 
represent the system adequately. 

Since many aspects of the landing-impact problem are so 
intimately connected with the mechanics of the landing gear, 
the subject of landing-gear behavior has received analytical 
treatment at various times (see bibliography). Many of 
the earlier investigations, in order to reduce the mathematical 
complexity of the analysis, were limited to consideration of 
highly simplified linear systems which have little relation to 
practical landing gears. Some of the more recent papers 
consider, with different degrees of simplification, more real- 
istic nonlinear systems. The present report represents an 
attempt at a more complete analysis of the mechanics of 
practical landing gears and, in addition, investigates the im- 
portance of the various elements which make up the landing 
gear, as well as the extent to which the system can be reason- 
ably simplified for the purpose of rapid analysis. 

The basic analysis is presented in a general form and takes 
into account such factors as the hydraulic (velocity square) 
resistance of the orifice, the forces due to air compression and 
internal friction in the shock strut, the nonlinear force- 
deflection characteristics of the tire, the wing lift, the inclina- 
tion of the landing gear and the effects of wheel spin-up drag 
loads. An evaluation of the applicability of the analysis to 
actual landing gears is presented for the case of a vertical 
landing gear in the absence of drag loads by comparing cal- 
culated results with drop-test data. 


Since some parameters, such as the dynamic force- 
deflection characteristics of the tire, the orifice discharge 
coefficient, and the polytropic exponent for the air-compression 
process, may not be accurately known in practical design 
problems, a study is made to assess the effects of variations 
in these parameters on the calculated landing-gear behavior. 

Studies are also presented to evaluate the extent to which 
the dynamical system can be simplified without greatly im- 
pairing the validity of the calculated results. In addition to 
the investigations for specific cases, generalized solutions for 
the behavior of a simplified system are presented for a wide 
range of landing-gear and impact parameters which may be 
useful in preliminary design. 

SYMBOLS 


A a 

Ah 

A 0 
A i 

A2 

A P 


C d 

d 

Fa 

Fn 

F r 

Fs 
F 1 

F 2 


Fy 

H 

F R 

Fs 

F n 



g 

Kl 

L 

U 


U 


a,b,m,r 

a' 

m' 

n 


pneumatic area 
hydraulic area 

area of opening in orifice plate 
internal cross-sectional area of shock-strut inner 
cylinder 

external cross-sectional area of shock-strut inner 
cylinder 

cross-sectional area of metering pin or rod in 
plane of orifice 
net orifice area 
orifice discharge coefficient 
overall diameter of tire 
pneumatic force in shock strut 
hydraulic force in shock strut 
friction force in shock strut 
total axial shock-strut force 
normal force on upper bearing (attached to inner 
cylinder) 

normal force on lower bearing (attached to outer 
cylinder) 

force normal to axis of shock strut, applied at 
axle 

vertical force, applied at axle 
horizontal force, applied at axle 
resultant force, applied at axle 
force parallel to axis of shock strut, applied to 
tire at ground 

force normal to axis of shock strut, applied to 
tire at ground 

vertical force, applied to tire at ground 
horizontal force, applied to tire at ground 
resultant force, applied to tire at ground 
gravitational constant 
lift factor, Lj\Y 
lift force 

axial distance between upper and lower bearings, 
for fully extended shock strut 
axial distance between axle and lower bearing 
(attached to outer cylinder), for fully ex- 
tended shock strut 

constants corresponding to the various regimes 
of the tire-deflection process 
combined constant, ad 
combined constant, md T 

polytropic exponent for air-compression process 
in shock strut 
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R 

Pa 

Ph 

0 

I'd 

S 

T 

t 

T 

V 

I w 

Vy 

V H 

W 

Wi 

w 2 

x 2 




z 2 


U 1 


U 2 

a 




Vs 


Reynolds number 

air pressure in upper chamber of shock strut 
hydraulic pressure in lower chamber of shock 
strut 

volumetric rate of discharge through orifice 

radius of deflected tire 

shock-strut axial stroke 

wheel inertia torque reaction 

time after contact 

time after beginning of shock-strut deflection 
air volume of shock strut 

polar moment of inertia for wheel assembly 
about axle 
vertical velocity 
horizontal velocity 
total dropping weight 
weight of upper mass above strut 
weight of lower mass below strut 
horizontal displacement of lower mass from 
position at initial contact 
vertical displacement of upper mass from posi- 
tion at initial contact 

vertical displacement of lower mass from posi- 
tion at initial contact 

dimensionless upper-mass displacement from 
position at initial contact 
dimensionless lower-mass displacement from 
position at initial contact 
dimensionless shock-strut stroke, u x -u 2 
dimensionless time after contact 
angle between shock-strut axis and vertical 

f u" da 

shock-strut effectiveness, 

'M'l mazO’max 


C U lmax , 

Jo Ul " dUl 

rij a landing-gear effectiveness, — ^ — 

no & & > u " u 

1 max l max 

e time interval in numerical integration procedures 

n coefficient of friction between tire and runway 

Hi coefficient of friction for upper bearing (attached 

to inner cylinder) 

ju 2 coefficient of friction for lower bearing (attached 

to outer cylinder) 

p mass density of hydraulic fluid 

a angular acceleration of wheel 

Axes: 

2 vertical axis, positive downward 

x horizontal axis, positive rearward 

Subscripts: 

0 at instant of initial contact 

r at instant of initial shock-strut deflection 

su at instant of wheel spin-up 

max maximum value 


Notation: 

|( )| absolute value of ( ) 

( )* estimated value of ( ) 

The use of dots over symbols indicates differentiation with 
respect to time t or r. 

Prime marks indicate differentiation with respect to 
dimensionless time 6. 


MECHANICS OF LANDING GEAR 
DYNAMICS OF SYSTEM 

In view of the fact that landing-gear performance appears 
to be relatively unaffected by the elastic deformations of 
the airplane structure (see, for example, refs. 1 and 2) par- 
ticularly since in many cases the main gears are located 
fairly close to the nodal points of the fundamental bending 
mode of the wing, that part of the airplane which acts on a 
given gear can generally be considered as a rigid mass. 
As a result, landing-gear drop tests are often conducted in 
a jig where the mass of the airplane is represented by a 
concentrated weight. In particular instances, however, such 
as in the case of airplanes having large concentrated masses 
disposed in an outboard position in the wings, especially 
airplanes equipped with bicycle landing gear, consideration 
of the interaction between the deformation of the airplane 
structure and the landing gear may be necessary to repre- 
sent the system adequately. 

Since the present report is concerned primarily with the 
mechanics of the landing gear, it is assumed in the analysis 
that the landing gear is attached to a rigid mass which has 
freedom only in vertical translation. The gear is assumed 
infinitely rigid in bending. The combination of airplane 
and landing gear considered therefore constitutes a system 
having two degrees of freedom (see fig. 1(a)) as defined by 
the vertical displacement of the upper mass and the vertical 
displacement of the lower or unsprung mass, which is also 
the tire deflection. The strut stroke s is determined by 
the difference between the displacements z x and z 2 and, in 
the case of inclined gears, by the angle <p between the axis 
of the strut and the vertical. For inclined gears, compression 
of the shock strut produces a horizontal displacement of the 
axle x 2 . From consideration of the kinematics of the system 

it can be seen that — — and x 2 =s sin <p= — z 2 )tan <p- 

COS ip 

In the analysis, external lift forces, corresponding to the 
aerodynamic lift, are assumed to act on the system through- 
out the impact. In addition to the vertical forces, arbitrary 
drag loads are considered to act between the tire and the 
ground. 

The system treated in the analysis may therefore be con- 
sidered to represent either a landing-gear drop test in a jig 
where wing lift and drag loads are simulated, or the landing 
impact of a rigid airplane if rotational motions are neglected. 
Rotational freedom of the airplane, where significant, may 
be taken into account approximately by use of an appro- 
priate effective mass in the analysis. 

Figure 1 (b) shows a schematic representation of a typical 
olco-pneumatic shock strut used in American practice. The 
lower chamber of the strut contains hydraulic fluid and the 
upper chamber contains air under pressure. The outer cyl- 
inder of the strut, which is attached to the upper mass, 
contains a perforated tube which supports a plate with a 
small orifice, through which the hydraulic fluid is forced to 
flow at high velocity as a result of the telescoping of the 
strut. The hydraulic pressure drop across the orifice thus 
produced resists the closure of the strut, and the turbulence 
created provides a powerful means of absorbing and dis- 
sipating a large part of the impact energy. In some struts 
the orifice area is constant; whereas, in other cases a metering 
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(b) 


(a) System with two degrees of freedom. 


(b) Schematic representation of shock strut. 


Figure 1 . — Dynamical system considered in analysis. 


pill or rod is used to control the size of the orifice and govern 
the performance of the strut. 

The compression of the strut produces an increase in the 
air pressure which also resists the closure of the strut. In 
figure 1 (b) p h represents the oil pressure in the lower chamber 
and p a represents the air pressure in the upper chamber. 

In addition to the hydraulic resistance and air-pressure 
foices, internal bearing friction also contributes forces which 
can appreciably affect the behavior of the strut. 

The forces created within the strut impart an acceleration 
to the upper mass and also produce an acceleration of the 
lower mass and a deflection of the tire. Figure 1 (c) shows 
the balance of forces and reactions for the wheel, the inner 
cylinder, and the outer cylinder. It is clear that the strut 


and the tire mutually influence the behavior of one another 
and must be considered simultaneously in analyzing the system. 

FORCES IN SHOCK STRUT 

From consideration of the pressures acting in the shock 
strut it can be readily seen from figure 1 (b) that the total 
axial force due to hydraulic resistance, air compression, and 
bearing friction can be expressed by 

Fs—Vh(^\~ + Pa {A2—A1 )-\'p a A p J r F f 

where 

A 1 internal cross-sectional area of inner cylinder 

A 2 external cross-sectional area of inner cylinder 

A v cross-sectional area of metering pin or rod in plane 
of orifice 
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(c) Balance of forces and reactions for landing-gear components. 
Figure 1 . — Concluded. 
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This expression can also be written as 


Fs— ijph — Pa) (A i — A p )-\-p a Az-\-F f 
— ( Vh~Va)A h -\-p a A a -\-F f 

= F *+ F a -\- F f ( 1 ) 

where 

Vh—Pa pressure drop across the orifice 

A h hydraulic area (A { —A p for the strut shown in fig. 1) 

A a pneumatic area ( A 2 for the strut shown in fig. 1) 

In this report the terms (p h —p a )A h and p a A a are referred 
to as hydraulic force F h and pneumatic force F a , respec- 
tively. For the strut shown in figure 1, the hydraulic 
and pneumatic areas are related to the strut dimensions as 
previously noted. In the case of struts having different 
internal configurations, the hydraulic and pneumatic areas 
may bear somewhat different relations to the dimensions 
of the strut. In such cases, however, consideration of the 
pressures acting on the various components of the strut 
should permit these areas to be readily defined. 

Hydraulic force. — The hydraulic resistance in the shock 
strut results from the pressure difference associated with the 
flow through the orifice. In a landing gear the orifice area 
is usually small enough in relation to the diameter of the 
strut so that the jet velocities and Reynolds numbers are 
sufficiently large that the flow is fully turbulent. As a 
result the damping force varies as the square of the tele- 
scoping velocity rather than linearly with the velocity. 
Since the hydraulic resistance is the major component of 
the total shock-strut force, viscous damping cannot be 
reasonably assumed, even though such an assumption 
would greatly simplify the analysis. 

The hydraulic resistance can be readily derived by making 
use of the well-known equation for the discharge through 
an orifice, namely, 

Q = C d A n ^l (p h —p a ) 

where 

0 volumetric rate of discharge 
C d coefficient of discharge 
A n net orifice area 

p h hydraulic pressure in lower chamber 
p a aii* pressure in upper chamber 
p mass density of hydraulic fluid 
From considerations of continuity, the volumetric rate of 
discharge can also be expressed as the product of the tele- 
scoping velocity s and the hydraulic area A h 

Q — A h s 


Equating the preceding expressions for the discharge per- 
mits writing the following simple equation for the pressure 
drop across the orifice 


Ph~Pa = 


pA h 2 s 2 
2 (C d A n y 


Ph—pa and the area A h which is subjected to the hydraulic 
pressure, as previously noted. Thus 


F h =- ( 


pA h 


'2 (C a A„y s * (2) 

Equation (2) can be made applicable to both the compres- 

s 

sion and elongation strokes by introducing the factor y-r-j 
to indicate the sign of the hydraulic resistance; thus 


V — 1 pAI. 
h \s\2(C d A n y 


(2a) 


The net orifice area A n may be either a constant or, when a 
metering pin is used, can vary with strut stroke; that is, 
A n =A 0 — A p =A n (s), where A 0 is the area of the opening in 
the orifice plate and A p is the area of the metering pin in the 
plane of the orifice. At the present time there appears to be 
some tendency to eliminate the metering pin and use a con- 
stant orifice area, particularly for large airplanes, in which 
case A n =A 0 . In the general case, the orifice discharge 
coefficient might be expected to vary somewhat during an 
impact because of changes in the size and eonfiguration of 
the net orifice area, changes in the exit conditions on the 
downstream face of the orifice due to variations in the amount 
of hydraulic fluid above the orifice plate, changes in the entry 
conditions due to variations in the length of the flow chamber 
upstream of the orifice, and because of variations in the 
Reynolds number of the flow, so that, in general, C d =C d (s , R). 
Although the individual effects of these factors on the dis- 
charge coefficients for orifices in shock struts have not been 
evaluated, there is some experimental evidence to indicate 
appreciable variations of the discharge coefficient during 
impact, particularly in the case of struts with metering pins. 
It might be expected that such variations would be con- 
siderably smaller for gears having a constant orifice area. 

In order to evaluate the precision with which the orifice 
discharge coefficient has to be known, a brief study is 
presented in a subsequent section which shows the effect of 
the discharge coefficient on the calculated behavipr of a 
landing gear with a constant orifice area, under the assump- 
tion that the discharge coefficient is constant during the 
impact. 

The foregoing discussion has been concerned primarily with 
the compression stroke of the shock strut. Most struts 
incorporate some form of pressure-operated rebound check 
valve, sometimes called a snubber valve, which comes into 
action after the maximum stroke has been attained and closes 
off the main orifice as soon as the strut begins to elongate, so 
that the fluid is forced to return to the lower chamber through 
small passages. The action of the snubber valve introduces 
greatly increased hydraulic resistance to dissipate the energy 
stored in the strut in the form of air pressure and to prevent 
excessive rebound. The product C d A n to be used in equation 
(2a) during the elongation stroke is generally uncertain. The 
exact area A n during elongation is usually somewhat difficult 
to define from the geometry of the strut since in many cases 
the number of connecting passages varies with stroke and the 


The hydraulic resistance F h due to the telescoping of the 
strut is given by the product of the differential pressure 
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leakage area around the piston may be of the same order of 
magnitude as the area of the return passages. Furthermore, 
the magnitude of the orifice discharge coefficient, and even 
possibly the nature of the resistance, are questionable due to 
the foaming state of the returning fluid. Fortunately, the 
primary interest is in the compression process rather than 
the elongation process since the maximum load always occurs 
before the maximum strut stroke is reached. 

Pneumatic force. — The air-pressure force in the upper 
chamber is determined by the initial strut inflation pres- 
sure, the area subjected to the air pressure (pneumatic area), 
and the instantaneous compression ratio in accordance 
with the poly tropic law for compression of gases, namely 
p a v n = Constant, or 



where 

p a air pressure in upper chamber of shock strut 

Pao air pressure in upper chamber for fully extended strut 

v air volume of shock strut 

v 0 air volume for fully extended strut 

Since the instantaneous air volume is equal to the difference 

between the initial air volume and the product of the stroke 


and pneumatic area A a , 



The force due 


to the air pressure is simply the product of the pressure and 
the pneumatic area: 


(3> 


In the preceding equations, the effective poly tropic 
exponent n depends on the rate of compression and the rate 
of heat transfer from the air to the surrounding environment . 
Low rates of compression would be expected to result in 
values of n approaching the isothermal value of 1.0; whereas 
higher values of n, limited by the adiabatic value of 1.4, 
would be expected for higher rates of compression. The 
actual thermodynamic process is complicated by the violent 
mixing of the highly turbulent efflux of hydraulic fluid and 
the air in the upper chamber during impact. On the one 
hand, the dissipation of energy in the production of turbu- 
lence generates heat; on the other hand, heat is absorbed by 
the aeration and vaporization of the fluid. The effect of this 
mixing phenomenon on the poly tropic exponent or on the 
equivalent air volume is not clear. A limited amount of 
experimental data obtained in drop tests (refs. 3 and 4), 
however, indicates that the effective poly tropic exponent 
may be in the neighborhood of 1.1 for pz-actical cases. A 
brief study of the importance of the air-compression process 
and the effects which different values of n may have on the 
calculated behavior of the landing gear is presented in a 
subsequent section. 

Internal friction force. -In the literature on machine de- 
sign the wide range of conditions under which frictional 
resistance can occur between sliding surfaces is generally 
classified in three major categories, namely, friction between 
dry surfaces, friction between imperfectly lubricated surfaces, 
and friction between perfectly lubricated surfaces. In the 


case of dry friction, the resistance depends on the physical 
characteristics of the sliding surfaces, is essentially propor- 
tional to the normal force, and is approximately independent 
of the surface area. The coefficient of friction /z, defined as 
the ratio of the frictional resistance to the normal force, is 
generally somewhat greater under conditions of rest (static 
friction) than under conditions of sliding (kinetic friction). 
Although the coefficient of kinetic friction generally de- 
creases slightly with increasing velocity, it is usually con- 
sidered, in first approximation, to be independent of velocity. 
If, on the other hand, the surfaces are completely separated 
by a fluid film of lubricant, perfect lubrication is said to exist. 
Lender these conditions the resistance to relative motion 
depends primarily on the magnitude of the relative velocity, 
the physical characteristics of the lubricant, the area, and 
the film thickness, and is essentially independent of the 
normal force and the characteristics of the sliding surfaces. 
Perfect lubrication is rarely found in practice but is most 
likely under conditions of high velocity and relatively small 
normal pressure, where the shape of the sliding surfaces is 
conducive to the generation of fluid pressure by hydro- 
dynamic action. In most practical applications involving 
lubrication, a state of imperfect lubrication exists and the 
resistance phenomenon is intermediate between that of dry 
friction and perfect lubrication. 

In the case of landing-gear shock struts, the conditions 
under which internal friction is of concern usually involve 
relatively high normal pressures and relatively small sliding 
velocities. Moreover, the usual types of hydraulic fluid 
used in shock struts have rather poor lubricating properties, 
and the shape of the bearing surfaces is generally not con- 
ducive to the generation of hydrodynamic pressures. It 
would therefore appear that the lubrication of shock strut 
bearings is, at best, imperfect; in fact, the conditions appear 
to approach closely those for dry friction. In the present 
analysis, therefore, it is assumed, in first approximation, that 
the internal friction between the bearings and the cylinder 
walls follows laws similar to those for dry friction; that is, 
the friction force is given by the product of the normal force 
and a suitably chosen coefficient of friction. 

With these assumptions the internal friction forces pro- 
duced in the strut depend on the magnitude of the forces on 
the axle, the inclination of the gear, the spacing of the bear- 
ings, and the coefficient of friction between the bearings and 
the cylinder walls. Figure 1(c) schematically illustrates the 
balance of forces acting on the various components of the 
landing gear. The total axial friction in the shock strut is 
the sum of the friction forces contributed by each of the 
bearings: 

where 

F f axial friction force 

Mi coefficient of friction for upper bearing (attached to 
inner cylinder) 

F x normal force on upper bearing (attached to inner 
cylinder) 

M 2 coefficient of friction for lower bearing (attached to 
outer cylinder) 
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F 2 normal force on lower bearing (attached to outer 
cylinder) 

t^t factor to indicate sign of friction force 

\ s \ 

During the interval prior to the beginning of shock-strut 
motion the friction forces depend on the coefficients of static 
friction; after the strut begins to telescope the coefficients of 
kinetic friction apply. 

From considerations of the balance of moments it can be 
seen from figure 1 (c) that 



so that 

[<*"+«> It5+«] «> 

where 

F N(t = F Va sin <p— F„ a cos <p (4a) 

and 

F Na force normal to strut applied at axle 
F Va vertical force applied at axle 
F Hq horizontal force applied at axle 

angle between strut axis and vertical 
/ j axial distance between upper and lower bearings, for 

fully extended strut 

/ 2 axial distance between axle and lower bearing (attached 
to outer cylinder), for fully extended strut 
The quantities F Na , F Va , and F Ha are forces applied at the 
axle and differ from the ground reactions by amounts equal 
to the inertia forces corresponding to the respective accelera- 
tion component of the lower mass. Since the inner cylinder 
generally represents only a relatively small fraction of the 
lower mass, the lower mass may reasonably be assumed to be 
concentrated at the axle. With this assumption, the rela- 
tionships between the forces at the axle and the forces at the 
ground are given by 



The normal force at the axle can therefore be expressed in 
terms of the ground reactions and the component accelera- 
tions of the lower mass by 

Fff a =(F Vg + ** 2 1 2 — W r 2 ) sin <p — ( F„— ** ? i 2 ) t'os ip (4b) 
where 

F v vertical force applied to tire at ground 
F Iffi horizontal force applied to tire at ground 
11 2 effective mass below shock strut, assumed concentrated 
g at axle 

x 2 horizontal acceleration of axle 

z 2 vertical acceleration of axle 

In the case of an inclined landing gear having infinite stiff- 
ness in bending, the horizontal displacement of the lower mass 
x 2 is related to the vertical displacements of the upperand lower 


masses by the kinematic relationship x 2 = (zi—z 2 ) tan <p, as 
previously noted. Double differentiation of this relation- 
ship gives (Si — 2 2 ) tan <p. Substitution of this expression 
into equation (4b) gives 


W 

F n —F v sin <p—F H cos — z x sin <p — W 2 sin <p (4c) 


In equation (4c) the quantity z { sin (p represents the ac- 
celeration of the lower mass normal to the strut axis when the 
gear is rigid in bending. In the case of a gear flexible in 
bending, the normal acceleration of the lower mass is not 
completely determined by the vertical acceleration of the 
upper mass and the angle of inclination of the gear. If it 
should be necessary to take into account, in particular cases, 
the effects of gear flexibility on the relationship between the 
normal force on the axle and the ground reactions, the quan- 
tity Z\ sin <p in equation (4c) may be replaced by estimated 
values of the actual normal acceleration of the lower mass as 
determined from consideration of the bending response of the 
gear to the applied forces normal to the gear axis. The effects 
of gear flexibility are not considered in more detail in the 
present analysis. 


FORCES ON TIRE 

Figure 2 (a) shows dynamic force-deflection characteris- 
tics for a 27-inch smooth-contour (type I) tire inflated to 32 
pounds per square inch. These characteristics were deter- 
mined from time-history measurements of vertical ground 
force and tire deflection in landing-gear drop tests with a 
nonrotating wheel at several vertical velocities. As can be 
seen, the tire compresses along one curve and unloads along 
another, the hysteresis loop indicating appreciable energy 
dissipation in the tire. There is some question as to whether 
the amount of hysteresis would be as great if the tire were 
rotating, as in a landing with forward speed. The force- 
deflection curve for a velocity of 11.63 feet per second is for 
a severe impact in which tire bottoming occurs and shows 
the sharp increase in force with deflection subsequent to 
bottoming. 

In figure 2 (b) the same force-deflection characteristics 
are shown plotted on logarithmic coordinates. As can be 
seen, the force exhibits an exponential variation with deflec- 
tion. A systematized representation of the force-deflection 
relationship can therefore be obtained bv means of simple 
equations having the form 



where 

F Vg vertical force, applied to tire at ground 
z 2 vertical displacement of lower mass from position at 
initial contact (radial deflection of tire) 
d overall diameter of tire 

m, r constants corresponding to the various regimes of the 
tire-deflection process 
m' combined constant, md r 
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(a) Uniform coordinates. (b) Logarithmic coordinates. 

Figure 2. — Dynamic force-deflection characteristics of tire. 


It may be noted from figure 2 that essentially the same 
force-deflection curve holds during compression for all impact 
velocities, up to the occurrence of tire bottoming, and that in 
figure 2 (b) the slopes of the curves in each of the several 
regimes of the tire-deflection process are also independent of 
velocity, except in the compression regime following tire 
bottoming. 

Figure 2 also shows simple approximations to the tire 
characteristics which were obtained by fitting straight-line 
segments (long-daslied lines) to the actual force-deflection 
curves in figure 2 (a) for impacts at 8.86 and 11.63 feet per 
second. These approximations, hereinafter referred to as 
linear-segment approximations, are included in a study, 
presented in a subsequent section, to evaluate the degree of 
accuracy required for adequate representation of the tire 
characteristics. The various representations of the tire 
characteristics considered and the pertinent constants for 
each regime of tire deflection are shown in figure 3. 

EQUATIONS OF MOTION 

The internal axial force F s produced by the shock strut 
was shown in a previous section to be equal to the sum of the 


hydraulic, pneumatic, and friction forces, as given by 
equation (1). Since these forces act along the axis of the 
strut, which may be inclined to the vertical by an angle <p, 
the vertical component of the axial shock-strut force is given 
by F s cos <p. The vertical component of the force normal 
to the shock strut is given by F Na sin <p. These forces act 
in conjunction with the lift force and weight to produce an 
acceleration of the upper mass. The equation of motion 
for the upper mass is 


F s cos (p-\-F Na sin p+L — W\= 



( 6 ) 


The vertical components of the axial and normal shock- 
strut forces also act, in conjunction with the weight of the 
lower mass, to produce a deformation of the tire and an 
acceleration of the lower mass. The equation of motion 
for the lower mass is 


F s cos <p~\~ Ftf a sin <p + W 2 — z 2 =F Vg (z 2 ) (7) 


where the vertical ground reaction F Vg is expressed as a 
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Tire deflection, z 2 , ft 

(a) Impact without tire bottoming, IV 0 =8.86 feet per second. (b) Impact with tire bottoming, F K() = 11.63 feet per second. 

Figure 3. — Tire characteristics considered in solutions (logarithmic coordinates). 


function of the tire deflection z 2 . The relationship between 
F v g and z 2 has been discussed in the previous section on tire 
characteristics. 

By combining equations (6) and (7), the vertical ground 
force can be written in terms of the inertia reactions of the 
upper and lower masses, the lift force, and the total weight. 
The overall dynamic equilibrium is given by 

FvM=-—h-^r*t-L+w ( 8 ) 

* g g 

MOTION PRIOR TO SHOCK-STRUT DEFLECTION 

Conventional oleo-pneumatic shock struts are inflated 
to some finite pressure in the fully extended position. Thus 
the strut does not begin to deflect in an impact until sufficient 
force is developed to overcome the initial preloading imposed 
by the air pressure and internal friction. Since the strut is 
effectively rigid in compression, as well as in bending, 
prior to this instant, the system may be considered to have 
only one degree of freedom during the initial stage of 


the impact. The equations of motion for the one-degree-of- 
freedom system are derived in order to permit determination 
of the initial conditions required for the analysis of the 
landing-gear behavior subsequent to the beginning of shock- 
strut deflection. 

Since z x =z 2 =z during this first phase of the impact, 
equation (8) may be written as 

FvAz)=~2-W(K l -1) (9) 

* g 

where 



For the general case of an exponential relationship between 
vertical ground force and tire deflection, equation (5) applies 
and the equation of motion becomes 

— S + m2 r 4- W(K L — 1) = 0 

g 


( 10 ) 
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The shock strut begins to telescope when the sum of the inertia, weight, and lift forces becomes equal to the vertical 
components of the axial and normal shock-strut forces. At this instant t T , F s =F ao +F, r and equation (6) can be written as 

( F„ o + F !r ) cos <p + F N ^ sin <p+K l H '—IT, 


F„ n initial air-pressure preload force, p„ 0 A a 
F fr static friction at instant t T 


At the instant t T , s = 0 and equation (4) becomes 


Fr = F Na 'K n 


where 


(11a) 




and fj . ! and ju 2 are coefficients of static friction. 

wSince the strut is assumed essentially rigid in compression (and also rigid in bending), there is no kinematic displ„„ 
ment of the lower mass in the horizontal direction up to the beginning of shock-strut deflection, so that x 2 =0 and equation 
(4b) becomes 


ace- 


FN a=(^ F Vg T +^Y ^ T ~ W z) sin V~ F H gT COS 

Incorporating equations (11a), (lib), and (9) into equation (11) gives 


(lib) 


.. ^ a 0 ( i Kn sin <p cos <p ) (Al W Wi) /* h i K ^ cos tp -|- sin ip) 

Zt ~W 

( ± sin ip — cos tp) 


( 12 ) 


In equation (12) wherever the ± sign appears, the plus signs 
apply when F Nq ^> 0 and the minus signs apply when F N(i <0. 

From equation (10) the vertical displacement of the system 
at the instant t T is given in terms of the corresponding 
acceleration by 

(,3) 

Integrating equation (10) and noting that z 0 = 0 provides 
the relationship between the vertical velocity and the 
vertical displacement of the system at the beginning of 
shock-strut deflection 

Zr=yj ^« 2 -f z/ +1 + W{K l -\)z^ (14) 


where the general expression for the variable z is obtained 
from equation (14) without the subscripts r. Performing 
the indicated integration gives 

tr= M% l sin_1 CCl-^-sin- 1 (15) 

where 

-i g = 

V*o*^+[(l ~K L )g]* 

The computation of t T can be greatly simplified by use of 
the following approximation which assumes a linear relation- 
ship between velocity and time: 


In view of the fact that the tire force-deflection curve is 
essentially linear for small deflections, it may be reasonably 
assumed that r= 1 for the purpose of determining the time 
after contact at which the strut begins to telescope. With 
this assumption t T can be determined from the relationship 



tr 


2Z T 


(15a) 


Equation (15a) should be a fairly good approximation in 
view of the relatively short time interval between initial con- 
tact and the beginning of shock-strut motion. 

Equations (12), (13), and (14) permit the determination 
of the vertical acceleration, displacement, and velocity, re- 
spectively, of the system (upper and lower masses) at the 
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beginning of shock-strut deflection. Equation (15) or (15a) permits calculation of the time interval between initial contact 
and this instant. These equations provide the initial conditions required for the analysis of the behavior of the landing 
gear as a system with two degrees of freedom after the shock strut begins to deflect. 

If drag loads are considered, the solution of equation (12) requires knowledge of the horizontal ground force F„^ at 
the instant t T . Since the present analysis does not explicitly treat the determination of drag loads, values of /' H ^ have 
to he estimated, either from other analytical considerations, experimental data, or on the basis of experience. 


MOTION SUBSEQUENT TO BEGINNING OF SHOCK-STRUT DEFLECTION 

Once the sum of the inertia, weight, and lilt forces becomes sufficiently large to overcome the preloading force in the 
shock strut due to initial air pressure and internal friction, the shock strut can deflect and the system becomes one having 
two degrees of freedom. Incorporating the expressions for the hydraulic, pneumatic, and friction forces (eqs. (2a), (3), 
and (4)) into equation (6) permits the equation of motion for the upper mass to be written as follows: 


Wi 


2i + 


{w\ 2 (OA.y bi+PaoAa (v 0 -A a s) + | S , l-.lt 


(Ml + Mi) 


l 2 — s 

l\-\-S 


/X2JI COS <p-\- KlW — WiH ^ — 0 


(16) 


where 

Z\ — Z 2 

$ = 

COS <p 


s = 


Z\ — Zo 
COS <p 


and, since F v * =F Vg (z 2 ), equation (4c) becomes 


F n —F v (z 2 ) sin <p— F h cos<p+ — -Z\ sin <p—W 2 sin <p 

a g g g 

where F v (z 2 ) is determined from the force-deflection characteristics of the tire. For the usual type of pneumatic tire, 
F v (z 2 )=mz 2 r , as previously noted. 

Similarly, the equation of motion for the lower mass follows from equation (7): 

cos <H- F v (z 2 )— Fn sin <p—W 2 = 0 (17) 


T Mpi (spfe)' + UFi +M ]} 

The overall dynamic equilibrium equation is still, of course, as given by equation (8) 


1 2 i+ " 22+11 (Kl — Fy g (z 2 ) — 0 


Any two of the preceding equations (eqs. (16), (17), and 
(8)) are sufficient to describe the behavior of the landing 
gear subsequent to the beginning of shock-strut motion. 
These equations may be used to calculate the behavior of a 
given landing-gear configuration or to develop orifice and 
metering-pin characteristics required to produce a specified 
behavior for given impact conditions. They may also be 
used as a basis for the calculation of dynamic loads in flexible 
airplane structures either by (a) determining the landing- 
gear forcing function under the assumption that the upper 
mass is a rigid body and then using this forcing function to 
calculate the response of the elastic system or (b) combining 
the preceding equations with the equations representing the 
additional degrees of freedom of the structure; the simul- 
taneous solution of the equations for such a system would 
then take into account the interaction between the deforma- 
tion of the structure and the landing gear. 


SOLUTION OF EQUATIONS OF MOTION 

In the general case the analysis of a landing gear involves 
the solution of the equations of motion given in the section 
entitled “Motion Subsequent to the Beginning of Shock- 
Strut Deflection,” with the initial conditions taken as the 
conditions of motion at the beginning of shock-strut deflec- 
tion, as determined in accordance with the initial impact 
conditions and the equations given in the section entitled 
“Motion Prior to Shock-Strut Deflection.” 

NUMERICAL INTEGRATION PROCEDURES 

In view of the fact that the equations of motion for the 
landing gear subsequent to the beginning of shock-strut 
deflection are highly nonlinear, analytical solution of these 
equations does not appear feasible. In the present report, 
therefore, finite-difference methods are resorted to for the 
step-by-step integration of the equations of motion. Al- 
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though such numerical methods lack the generality of ana- 
lytical solutions and are especially time consuming if the 
calculations are carried out manually, the increasing availa- 
bility of automatic calculating machines largely overcomes 
these objections. 

Most of the solutions presented in this report were obtained 
with a procedure, hereinafter referred to as the “linear pro- 
cedure/ 7 which assumes changes in the motion variables to 
be linear over finite time intervals. A few of the solutions 
presented were obtained with a procedure, hereinafter referred 
to as the “quadratic procedure/ 7 which assumes a quadratic 
variation of displacement with time for successive intervals. 
The generalized solutions for the simplified equations dis- 
cussed in a subsequent section were obtained by means of 
the Runge-Kutta procedure. The application of these 
procedures is described in detail in appendix A. 

USE OF TIRE FORCE-DEFLECTION CHARACTERISTICS 

In order to obtain solutions for particular cases, it is, of 
course, necessary to have, in addition to information regard- 
ing the physical characteristics of the landing gear, some 
knowledge of the force-deflection characteristics of the tire. 

If extensive data regarding the dynamic tire character- 
istics, such as shown in figures 2 and 3, are available, an 
accurate solution can be obtained which takes into account 
the various breaks in the force-deflection curves (logarithmic 
coordinates), as well as the effects of hysteresis. In view of 
the fact that the constants m' and r have the same values 
throughout practically the entire tire compression process 
regardless of the initial impact velocity or the maximum 
load attained, these values of m r and r, as determined from 
the force-deflection curves, can be used in the calculation of 
the motion subsequent to the beginning of shock-strut deflec- 
tion until the first break in the force-deflection curve is 
reached prior to the attainment of the maximum force. If 
the conditions for the calculations are the same as those for 
which force-deflection curves are available, the values of 
m' and r for each of the several regimes subsequent to the 
first break can also be determined directly from the force- 
deflection curves. In general, however, the conditions will 
not be the same and interpolation will be necessary to 
estimate the values of m' for the subsequent regimes. 
Such interpolation is facilitated, particular^ after the maxi- 
mum force-deflection point has been calculated, by the fact 
that each subsequent regime has a fixed value of r , regardless 
of the initial impact conditions. 

The use of the tire-deflection characteristics in the calcula- 
tions is greatly simplified if hysteresis is neglected since the 
values of m' and r which apply prior to the first break in the 
force-deflection curves are then used throughout the entire 
calculation, except in the case of severe impacts where tire 
bottoming occurs, in which case new values of m' and r are 
employed in the tire-bottoming regime. A similar situation 
exists with respect to the constants a' and b when the linear 
approximations which neglect hysteresis are used. These 
simplifications would normally be employed when only the 

287846—54 3 


tire manufacturer’s static or so-called impact load-deflection 
data are available, as is usually the case. 

EFFECT OF DRAG LOADS 

Although the present analysis permits taking into account 
the effects of wheel spin-up drag loads on the behavior of the 
landing gear, the determination of the drag-load time history 
is not treated explicitly. Thus, if it is desired to consider the 
effects of the drag load on the gear behavior, such as in the 
case of a drop test in which drag loads are simulated by 
reverse wheel rotation or in a landing with forward speed, it 
is necessary to estimate the drag load, either by means of 
other analytical considerations or by recourse to experimental 
data. As a first approximation the instantaneous drag force 
may be assumed to be equal to the vertical ground reaction 
multiplied by a suitable coefficient of friction /x ; that is, 
F Hg =F Vg n, up to the instant when the wheel stops skidding, 
after which the drag force may be assumed equal to zero. 
(The current ground-loads requirements specify a skidding 
coefficient of friction /x = 0.55; limited experimental evidence, 
on the other hand, indicates that ju may be as high as 0.7 or 
as low as 0.4.) In some cases experimental data indicate 
that representation of the drag-load time history can be 
simplified even further by assuming a linear variation of the 
drag force with time during the period of wheel skidding. 

The instant at which the wheel stops skidding can be 
estimated from the simple impulse-momentum relationship 




F v dt- 


IwVhq 
r d 2 


where 

I w polar moment of inertia of wheel assembly about axle 
V Hq initial horizontal velocity 
r d radius of deflected tire 
t su time of wheel spin-up 

When the drag force is expressed in terms of the vertical 

force, the value of the integral f F H dt can be determined as 

Jo g 

the step-by-step calculations proceed and the drag-force 
term eliminated from the equations of motion after the re- 
quired value of the integral at the instant of spin-up is 
reached. 

EVALUATION OF ANALYSIS BY COMPARISON OF CALCULATED 
RESULTS WITH EXPERIMENTAL DATA 

In order to evaluate the applicability of the foregoing 
analytical treatment to actual landing gears, tests were 
conducted in the Langley impact basin with a conventional 
oleo-pneumatic landing gear originally designed for a small 
military training airplane. A description of the test specimen 
and apparatus used is given in appendix B. 

In this section calculated results are compared with ex- 
perimental data for a normal impact and a severe impact 
with tire bottoming. The vertical velocities at the instant 
of ground contact used in the calculations correspond to 
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the vertical velocities measured in the tests. Equations 
(12), (13), (14), and (15a) were used to calculate the values 
of the variables at the instant of initial shock-strut deflection. 
Numerical integration of equations (16) and (17) provided 
the calculated results for the two-degree-of-freedom system 
subsequent to the beginning of shock-strut deflection. 

In these calculations the discharge coefficient for the orifice 
and the polytropic exponent for the air-compression process 
were assumed to have constant values throughout the impact. 
Consideration of the shape of the orifice and examination 
of data for rounded approach orifices in pipes suggested 
a value of C d equal to 0.9. Evaluation of data for other 
landing gears indicated that the air-compression process 
could be represented fairly well by use of an average value 
of the effective polytropic exponent n=1.12. In view of 
the fact that the landing gear was mounted in a vertical 
position and drag loads were absent in the tests, friction 
forces in the shock strut were assumed to be negligible 
in the calculations. Since the weight was fully balanced 
by lift forces in the tests, the lift factor K L was taken equal to 
1 . 0 . The appropriate exact tire characteristics (see fig. 3) 
were used for each case. 

NORMAL IMPACT 

Figure 4 presents a comparison of calculated results with 
experimental data for an impact without tire bottoming at 
a vertical velocity of 8.86 feet per second at the instant of 
ground contact. The exact dynamic force-deflection charac- 
teristics of the tire, including hysteresis, were used in the 
calculations. These tire characteristics are shown by the 
solid lines in figure 2 (a) and values for the tire constants 
in' and r are given in figure 3 (a). 

Calculated time histories of the total force on the upper 
mass and the acceleration of the lower mass are compared 
with experimental data in figure 4 (a). Similar comparisons 
for the upper-mass displacement, upper-mass velocity, lower- 
mass displacement, strut stroke, and strut telescoping 
velocity are presented in figure 4 (b). As can be seen, the 
agreement between the calculated and experimental results 
is reasonably good throughout most of the time history. 
Some of the minor discrepancies during the later stages of 
the impact appear to be due to errors in measurement since 
the deviations between the calculated and experimental 
upper-mass accelerations (as represented by the force on 
the upper mass) are incompatible with those for the upper- 
mass displacements, whereas the calculated upper-mass dis- 
placements are necessarily directly compatible with the 
calculated upper-mass accelerations. The maximum value 
of the experimental acceleration of the lower mass may be 
somewhat high because of overshoot of the accelerometer. 

In addition to the total force on the upper mass, figure 
4 (a) presents calculated time histories of the hydraulic 
and pneumatic components of the shock-strut force, as 
determined from equations (2) and (3), respectively. It 
can be seen that throughout most of the impact the force 
developed in the shock strut arises primarily from the hy- 
draulic resistance of the orifice. Toward the end of the 
impact, however, because of the decreased telescoping 
velocities and fairly large strokes which correspond to high 


compression ratios, the air-pressure force becomes larger 
than the hydraulic force. 

IMPACT WITH TIRE BOTTOMING 

Figure 5 presents a comparison of calculated and experi- 
mental results for a severe impact (V Vq = 11.63 ft per sec) 
in which tire bottoming occurred. The tire force-deflection 
characteristics used in the calculations are shown by the 
solid lines in figure 3 (b). Region (1) of the tire force- 
deflection curve has the same values of the tire constants 
m f and v as for the case previously discussed. Following 
the occurrence of tire bottoming, however, different values 
of m' and r apply. These values are given in figure 3 (b). 

It can be seen from figure 5 that the agreement between 
the calculated and experimental results for this case is 
similar to that for the comparison previously presented. 

The calculated instant of tire bottoming is indicated in 
figure 5. When tire bottoming occurs, the greatly increased 
stiffness of the tire causes a marked increase in the shock- 
strut telescoping velocity, as is shown in the right-hand 
portion of figure 5 (b). Since the strut is suddenly forced to 
absorb energy at a much higher rate, an abrupt increase 
in the hydraulic resistance takes place. The further increase 
in shock-strut force immediately following the occurrence of 
tire bottoming is evident from the left-hand portion of 
figure 5 (a). The sudden increase in lower-mass acceleration 
at the instant of tire bottoming can also be seen. 

In this severe impact the hydraulic resistance of the orifice 
represents an even greater proportion of the total shock- 
strut force than was indicated by the calculated results for 
an initial vertical velocity of 8.86 feet per second previously 
discussed. 

The foregoing comparisons indicate that the analytical 
treatment presented, in conjunction with reasonably straight- 
forward assumptions regarding the parameters involved in 
the equations, provides a fairly accurate representation of 
the behavior of a conventional oleo-pneumatic landing gear. 

PARAMETER STUDIES 

In the previous section comparisons of calculated results 
with experimental data showed that the equations which 
have been developed provide a fairly good representation of 
the behavior of the landing gear for the impact conditions 
considered. In view of the fact that the equations are 
somewhat complicated and require numerical values for 
several parameters such as the tire force-deflection constants 
m' and r, the orifice discharge coefficient C d , and the poly- 
tropic exponent n ) which may not be readily or accurately 
known in the case of practical engineering problems, it 
appears desirable (a) to determine the relative accuracy 
with which these various parameters have to be known and 
(b) to investigate the extent to which the equations can be 
simplified and still yield useful results. In order to accom- 
plish these objectives, calculations have been made to 
evaluate the effect of simplifying the force-deflection charac- 
teristics of the tire, as well as to determine the effects which 
different values of the orifice discharge coefficient and the 
effective polytropic exponent have on the calculated behavior. 
The results of these calculations are discussed in the present 
section. The question of simplification of the equations of 
motion is considered in more detail in a subsequent section. 


Velocity, fps 
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(a) Time histories of forces on upper mass and lower- mass acceleration. 

Figure 4. — Comparisons between calculated results and experimental data for normal impact; 
solution with exact exponential tire characteristics. VV 0 =8.86 feet per second; Cd = 0.9; ft = 1.12. 



(b) Time histories of landing-gear velocities and displacements. 
Figure 4. — Concluded. 
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(a) Time histories of forces on upper mass and lower-mass acceleration. 

Figure 5. — Comparisons between calculated results and experimental data for impact with tire bottoming; solution with 
exact exponential tire characteristics. F7 0 =11.63 feet per second; (7^=0. 9; w=1.12. 



(b) Time histories of landing-gear velocities and displacements. 
Figure 5. — Concluded. 
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REPRESENTATION OF TIRE FORCE-DEFLECTION CHARACTERISTICS 

In order to evaluate the degree of accuracy required for 
adequate representation of the tire force-deflection charac- 
teristics, comparisons are made of the calculated behavior 
of the landing gear for normal impacts and impacts with 
tire bottoming when the tire characteristics are represented 
in various ways. First, the force-deflection characteristics 
will be assumed to be exactly as shown by the solid-line 
curves in figure 2 (b), including the various breaks in the 
curve and the effects of hysteresis. These characteristics 
are referred to hereinafter as the exact exponential tire 
characteristics. The effects of simplifying the representa- 
tion of the tire characteristics will then be investigated by 
considering (a) the exponential characteristics without 
hysteresis; that is, the tire will be assumed to deflect and 
unload along the same exponential curve, (b) the linear- 
segment approximations to the tire characteristics (long- 
dashed lines), which also neglect hysteresis, and (c) errors 
introduced by neglecting the effects of tire bottoming in the 
case of severe impacts. The calculated results presented in 
this study make use of the relationships between vertical 
force on the tire and tire deflection, as shown in figures 
3 (a) and 3 (b). 

Figure 6 presents a comparison of the calculated results 
for a normal impact at a vertical velocity of 8.86 feet per 
second, whereas figure 7 permits comparison of the solutions 


for a severe impact, involving tire bottoming, at a vertical 
velocity of 11.63 feet per second. In figures 6 and 7 the 
solid-line curves represent solutions of the landing-gear 
equations when the exact exponential relationships between 
force and tire deflection are considered. Since these solu- 
tions were previously shown to be in fairly good agreement 
with experimental data (figs. 4 and 5), they are used as a 
basis for evaluating the results obtained when tire hysteresis 
is neglected and the force-deflection characteristics are repre- 
sented by either simplified exponential or linear-segment 
relationships. 

As in the calculations previously described, the solutions 
were obtained in two parts. During the first stage of the 
impact the shock strut was considered to be rigid until 
sufficient force was developed to overcome the initial air- 
pressure force. The calculations for the landing-gear behav- 
ior subsequent to this instant were based on the equations 
which consider the gear to have two degrees of freedom. 
Time histories of the upper-mass acceleration calculated on 
the basis of a rigid shock strut are shown by the dotted 
curves in figures 6 and 7. These solutions show the greatest 
rate of increase of upper-mass acceleration possible with 
the exponential tire force-deflection characteristics con- 
sidered. Comparison of these solutions with those for the 
two-degree-of-freedom system indicates the effect of the 
shock strut in attenuating the severity of the impact. 



(a) Time histories of upper-mass acceleration and lower-mass acceleration. 

Figure 6. Effect of tire characteristics on calculated landing-gear behavior in normal impact. FV 0 =8.86 feet per second; Cd— 0.9; n= 1.12, 
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(b) Time histories of landing-gear displacements and velocities. 
Figure 6. — Continued. 



(c) Time histories of shock-strut stroke and velocity, 
Figure 6. — Concluded. 


Displacement, ft w Upper- mass acceleration factor, >j /g 
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(a) Time histories of upper-mass acceleration and lower-mass acceleration. 

3 7. — Effect of tire characteristics on calculated landing-gear behavior for impact with tire bottoming. 
W 0 = 11.63 feet per second; Cd = 0.9; n = 1.12. 



(b) Time histories of landing^gear displacements and velocities. 
Figure 7. — Continued. 
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(c) Time histories of shock-strut stroke and velocity. 
Figure 7. — Concluded. 


Normal impact. — In the case of the normal impact at a 
vertical velocity of 8.86 feet per second, figure 6 shows that 
the solution obtained with the exponential force-deflection 
variation which neglects hysteresis and the solution with the 
linear-segment approximation to the tire characteristics are 
in fairly good agreement with the results of the calculation 
based on the exponential representation of the exact tire 
characteristics. The greatest differences between the solu- 
tions are evident in the time histories of upper-mass and 
lower-mass accelerations; considerably smaller differences are 
obtained for the lower-order derivatives, as might be ex- 
pected. With regard to the upper-mass acceleration, the 
three solutions are in very good agreement during the early 
stages of the impact. In the case of the simplified exponen- 
tial characteristics, neglect of the decreased slope of the 
force-deflection curve between the first break and the maxi- 
mum (regime 0 in fig. 3 (a)) resulted in the calculation of a 
somewhat higher value of the maximum upper-mass acceler- 
ation than was obtained with the exact tire characteristics. 
For the simplified exponential and linear-segment character- 
istics, neglect of hysteresis resulted in the calculation of 
somewhat excessive values of upper-mass acceleration sub- 
sequent to the attainment of the maximum vertical load. 
It is of interest to note that the calculated results for the 
exponential and linear-segment characteristics without hys- 
teresis were generally in quite good agreement with each 
other throughout the entire duration of the impact, although 
the assumption of linear-segment tire force-deflection char- 
acteristics did result in somewhat excessive values for the 
maximum lower-mass acceleration. On the whole, the 
simplified tire force-deflection characteristics considered per- 
mit calculated results to be obtained which represent the 
behavior of the landing gear in normal impacts fairly well. 


Impact with tire bottoming. — In the case of the severe 
impact at a vertical velocity of 11.63 feet per second, the 
effects of tire bottoming on the upper-mass acceleration, the 
lower-mass acceleration, and the strut telescoping velocity 
are clearly indicated in figure 7 by the calculated results 
based on the exact tire characteristics. As can be seen, the 
linear-segment approximation to the tire deflection character- 
istics which takes into account the effects of tire bottoming 
resulted in a reasonably good representation of the landing- 
gear behavior throughout most of the time history. On the 
other hand, as might be expected, the calculations which 
neglected the effects of bottoming on the tire force-deflection 
characteristics did not reveal the marked increase in the 
upper-mass acceleration due to the increased stiffness of the 
tire subsequent to the occurrence of bottoming. It is also 
noted that the discrepancies in the calculated upper-mass 
acceleration due to neglect of hysteresis in the later stages 
of the impact are more pronounced in this case than in the 
impact without tire bottoming previously considered, again 
as might be expected. 

EFFECT OF ORIFICE DISCHARGE COEFFICIENT 

In view of the fact that there is very little information 
available regarding the magnitude of discharge coefficients 
for orifices in landing gears, it appears desirable to evaluate 
the effect which differences in the magnitude of the orifice 
coefficient can have on the calculated results. Figure 8 
presents comparisons of calculated results for a range of 
values of the orifice discharge coefficient C d between 1.0 
and 0.7. The four solutions presented are for the same set 
of initial conditions as the normal impact without tire 
bottoming previously considered and are based on the 
exponential tire force-deflection characteristics which neglect 
hysteresis. 
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(a) Time histories of upper-mass acceleration and lower-mass acceleration. 


Figure 8. — Effect of orifice discharge coefficient on landing-gear behavior; calculations with exponential tire characteristics 

without hysteresis. Vr 0 =8.86 feet per second; w=1.12. 



(b) Time histories of landing-gear displacements and velocities. 
Figure 8. — Continued. 
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(c) Time histories of shock-strut stroke and velocity. 
Figure 8. — Concluded. 


Those calculations show that a decrease in the orifice dis- 
charge coefficient results in an approximately proportional 
increase in the upper-mass acceleration. This vari- 
ation is to be expected since the smaller coefficients cor- 
respond to reduced effective orifice areas which result in 
greater shock-strut forces due to increased hydraulic 
resistance. As a result of the increased shock-strut force 
acting downward on the lower mass, the maximum upward 
acceleration of the lower mass is reduced with decreasing 
values of the discharge coefficient. The increase in shock- 
strut force with decreasing discharge coefficient also results 
in a decrease in the strut stroke and telescoping velocity but 
an increase in the lower-mass velocity and displacement, 
as might be expected. However, since the increases in 
lower-mass displacement and velocity are smaller than the 
decreases in strut stroke and telescoping velocity, the upper- 
mass displacement and velocity are reduced with decreasing 
orifice discharge coefficient. 

These comparisons show that the magnitude of the orifice 
coefficient has an important effect on the behavior of the 
landing gear and indicates that a fairly accurate determi- 
nation of the numerical value of this parameter is necessary 
to obtain good results. 

EFFECT OF AIR-COMPRESSION PROCESS 

Since the nature of the air-compression process in a shock 
strut is not well-defined and different investigators have 
assumed values for the poly tropic exponent ranging any- 
where between the extremes of 1.4 (adiabatic) and 1.0 


(isothermal), it appeared desirable to evaluate the im- 
portance of the air-compression process and to determine the 
extent to which different values of the polytropic exponent 
can influence the calculated results. Consequently, solutions 
have been obtained for three different values of the poly- 
tropic exponent; namely, ti= 1.3, 1.12, and 0. 

The value n= 1.3 corresponds to a very rapid compression 
in which an adiabatic process is almost attained. The 
value n=1.12 corresponds to a relatively slow compression 
in which the process is virtually isothermal. The value 
r<, = 0 is completely fictitious since it implies constant air 
pressure within the strut throughout the impact. The 
assumption n = 0 has been considered since it makes one of 
the terms in the equations of motion a constant and permits 
simplification of the calculations. The three solutions 
presented are for the same set of initial conditions as the 
normal impact without tire bottoming previously con- 
sidered and are based on the exponential tire force-deflection 
characteristics which neglect hysteresis. 

Figure 9 shows that the air pressure contributes only a 
relatively small portion of the total shock-strut force through- 
out most of the impact since the compression ratio is rela- 
tively small until the later stages of the impact. Toward 
the end of the impact, however, the air-pressure force 
becomes a large part of the total force since the compression 
ratio becomes large, whereas the hydraulic resistance de- 
creases rapidly as the strut telescoping velocity is reduced 
to zero. 


Displacement, ft d Force on upper mass, lb 

f«v /n >■* • . Jo ' 
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(a) Time histories of upper-mass force and lower-mass acceleration. 

e 9. — Effect of polytropic exponent; calculations with exponential tire characteristics without hysteresis. 

W 0 = 8.86 feet per second; C d =0.9. 



(b) Time histories of landing-gear displacements and velocities. 
Figure 9. — Continued. 
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(c) Time histories of shock-strut stroke and velocity. 
Figure 9.- — Concluded. 


As a result, the calculations show that the magnitude of 
the poly tropic exponent has only a very small effect on the 
behavior of the landing gear throughout most of the impact. 
For the practical range of polvtropic exponents, variations 
in the air-compression process result in only minor differ- 
ences in landing-gear behavior, even during the very latest 
stages of the impact. The assumption of constant air 
pressure in the strut throughout the impact (n=0), however, 
does lead to the calculation of excessive values of stroke and 
of the time to reach the maximum stroke. The time history 
of the shock-strut force calculated on the basis of this 
assumption is, on the other hand, in quite good agreement 
with the results for the practical range of air-compression 
processes. 

On the whole it appears that the behavior of the landing 
gear is relatively insensitive to variations in the air- 
compression process. The foregoing results suggest that, in 
many cases, fairly reasonable approximations for the landing- 
gear force-time variation might be obtained even if the air- 
pressure term in the equations of motion were completely 
neglected. 

SIMPLIFICATION OF EQUATIONS OF MOTION 

The preceding studies have indicated that variations in the 
tire force-deflection characteristics and in the air-compression 
process individually have only a relatively minor effect on the 
calculated behavior of the landing gear. These results sug- 
gest that the equations of motion for the landing gear might 
be simplified by completely neglecting the internal air- 
pressure forces in the shock strut and by considering the tire 
force-deflection characteristics to be linear. With these as- 
sumptions, the equations of motion for the upper mass, 
lower mass, and complete system (eqs. (16), (17), and (8)) can 


be written as follows for the case where the wing lift is equal 
to the weight and the internal friction is neglected: 


— Zi + A(zi— z 2 ) 2 +W 2 =0 
g 


W 2 


Z 2 — A(z 1 — Z2) 2 -\-CbZ2-\~ b — 11 2 — 0 K 


Wi .. . w 2 .. . . 

— - z id z 2 -j- az 2 -\- b — 0 

s s 


(18) 


where 

, s pA h 3 
\s\ 2 (CdA n ) 2 cos<p 

and 

a slope of linear approximation to tire force-deflection 
characteristics 

b value of force corresponding to zero tire deflection, as 
determined from the linear-segment approximation to 
the tire force-deflection characteristics 
The motion variables at the beginning of shock-strut de- 
flection can be readily determined in a manner similar to 
that employed in the more general treatment previously dis- 
cussed. For the simplified equations the variables at the 
instant t T are given by 


z T = 



Z 




> 


(19) 


Z T = 


V 


z ° 


2 / 


Iii most cases the term -j|r z 2 is small in comparison with i 0 2 
so that z T ~z 0 . 
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(a) Time histories of upper-mass acceleration and lower-mass acceleration. 
Figure 10. — Evaluation of calculated results for simplified systems. W 0 = 8.86 feet per second; Cd=0.9. 


The values determined from equations (19) are used as 
initial conditions in the solution of equations (18). 

The fact that the lower mass is a relatively small fraction of 
the total mass suggests that the system might be simplified 
even further without greatly modifying the calculated results 
by assuming the lower mass to be equal to zero. With this 
assumption t T = 0 and the initial values of the variables in 
equations (18) correspond to the conditions at initial contact. 

EVALUATION OF SIMPLIFICATION 

In order to evaluate the applicability of these simplifica- 
tions, the behavior of the landing gear has been calculated in 
accordance with equations (18) for an impact with an initial 
vertical velocity of 8.86 feet per second. A similar calcula- 
tion has been made with the assumption I1 7 2 =0. These 
results are compared in figure 10 with the more exact solu- 
tions previously presented in figure 4, which include con- 
sideration of the air-compression springing and the exact 
exponential tire characteristics. A time history of the lower- 
mass acceleration is not presented for the case where W 2 is 
assumed equal to zero since the values of z 2 /g have no 
significance in this case. 

Figure 10 shows that the two simplified solutions are in 
quite good agreement with each other, as might be expected, 
and are also in fairly good agreement with the more exact 
results. Neglecting the air-pressure forces and assuming a 
linear tire force-deflection variation resulted in the calcula- 
tion of slightly lower values for the maximum upper-mass 
acceleration and somewhat higher values for the maximum 
stroke than were obtained with the more exact equations. 


The effect of neglecting the lower mass was primarily to 
reduce the lower-mass displacement (tire deflection), as a 
result of the elimination of the lower-mass inertia reaction. 

On the whole, it appears that the assumptions considered 
permit appreciable simplification of the equations of motion 
without greatly impairing the validity of the calculated results. 

GENERALIZED TREATMENT 

Equations and solutions. — By writing the simplified equa- 
tions of motion in terms of dimensionless variables, general- 
ized solutions can be obtained for a wide range of landing-gear 
and impact parameters which may be useful in pre- 
liminary design. If W 2 is taken equal to zero and it is 
further assumed that the tire force-deflection curve is 
represented by a single straight line through the origin 
(6=0 throughout the impact), equations (18) reduce to 

Zi + A(ii— z 2 ) 2 =0 j 

A(ii— z 2 ) 2 —az 2 =0 > (20) 

— zi+az 2 =0 \ 

g ) 

Wi 

where — > A } and a are constants, as previously defined, and 

any two of the foregoing equations are sufficient to describe 
completely the behavior of the system. With this representa- 
tion of the system, the shock strut begins to deflect at the in- 
stant of initial contact (t T = 0). Thus, the initial conditions 
for equations (20) are the initial impact conditions; namely, 

z 'o =Z2 o =0 and ^io =i 2 0 =2<,. 
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(b) Time histories of landing-gear displacements and velocities. 
Figure 10. — Continued. 



(c) 

j i 

.20 


(c) Time histories of shock-strut stroke and velocity. 
Figure 10. — Concluded. 
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As can be seen, with the equations in this form, the solution 

Wi 

depends on five parameters, namely, - — ■ > A, a, and the 


initial conditions z 0 and z 0 . However, since z 0 =0 in all 
cases, the number of variable parameters is reduced to four. 
In view of the fact that these parameters are independent of 
one another and each may take on a large range of values, a 
great many solutions and a large number of graphs would be 
required to cover the entire range of landing-gear and impact 
parameters with the equations in the form of equations (20). 

The number of independent parameters which have to be 
considered may be greatly decreased by the introduction of 
generalized dimensionless variables and the corresponding 
transformations of equations (20). In this case, generalized 
variables can be obtained which permit transformation of the 
equations of motion to a form which does not involve any 
constants. With the equations in this form, there is only 
one variable parameter, namely, the initial velocity param- 
eter. To determine the generalized variables which satisfy 
the aforementioned requirements, let 



u=za. 

and 

o=tp 

Thus, 

II 

II 

> 

and 

,, du’ 

u ~ de~ 


Substituting these new variables permits equations (20) 
to be written as 


(u/— « 2 ') 2 +(^ ^) Wi"=0 

(w/-- a 2 ') 2 -(^ jp) «2=0 

Ul " + (wji) 112=0 


(20a) 


The number of independent parameters will be reduced 
if all the combined constants in equations (20a) are set 
equal to one another, that is, let 

a W\_ a a _ a/ A 2 

aY~aJ 2 ~W^ 


From this relationship, it can be seen that 

A 


and 

a= 

Wt/g 


0=1 

Jwjg 

Thus the generalized variables become 

“■ =zi G 4 ) 


$ 

II 

, dui 

Ui =w =2 n 

lA 2 /a 

VWr/g 

, du 2 
W2= W 


r#Ta 

Wt/g 


and 



where 

With these new variables equations (20) can be written as 2 

(ui —U2') 2 +Ui ff =0\ 

(u 1 '—u 2 'y—u 2 =0 > (21) 

Ui"+u 2 =0 ) 

where any two of these equations are sufficient to describe 
the behavior of the system. 

Inasmuch as equations (21) do not involve any constants, 
their solutions are completely determined by the initial val- 
ues of the variables. Since the displacements at initial 
contact Ui Q and u 2q are equal to zero and the initial velocities 
Ui Q ' and u 2q ' are equal, the only parameter is the initial 
dimensionless velocity 


U Q '=Zo 


I A 2 g 
V Wia 


where u 0 ' =Ui 0 ' =u 2q ' . 

Generalized solutions of equations (21) are presented in 
figure 11 for values of Uq corresponding to a wide range of 
landing-gear and impact parameters. Parts (a) to (e) of 
figure 11 show the variations of the dimensionless variables 
during the impact; parts (f) and (g) show the maximum 
values of the more important variables as functions of Uq . 
Part (h) shows the shock-strut effectiveness tj s and the 
landing-gear effectiveness rji g . The shock-strut effective- 
ness, sometimes called “efficiency” and, in Europe, “plani- 
metric ratio,” is defined as 

C^maz 
J Ui" da 

U " max Cmax 

where a=U\— u 2 is the dimensionless shock-strut stroke. 
Since rj s represents the ratio of the energy actually absorbed 
by the shock strut to the maximum energy which the strut 
could possibly absorb for any combination of maximum 
acceleration (or load) and maximum stroke, it serves as a 
measure of the extent to which a given combination of 
maximum load and stroke has been utilized to absorb the 
energy of an impact. A similar measure of the energy 
absorption effectiveness of the landing gear as a whole is 
given by r) lg which is defined by 


Vig= 


L 


U lmaz 

U\" dui 


Ui U\ 

1 max l max 


2 Equations (21) maybe reduced to a single equation in one variable by differentiating the 
last equation and substituting for m' in the first equation. This gives 
By introducing the new variable w=u/ t this equation may be reduced to the second- 
order equation (w+u>")H-w'=0, subject to the initial conditions wo=uo' and uY= 0. 


Dimensionless upper-moss displocement, 
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Dimensionless time, 9 

(a) Relationship between upper-mass acceleration, lower-mass displacement, and time. 

Figure 11. — Generalized solutions for simplified system. 



(b) Relationship between upper-mass displacement and time. 
Figure 11. — Continued. 
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(c) Relationship between upper-mass velocity, lower-mass velocity, and time. 
Figure 11. — Continued. 


The generalized results presented in figure 11 3 can be 
used to estimate the performance of a given landing gear of 
known configuration for particular impact conditions or to 
choose the dimensions for a landing gear when the impact 
conditions and desired performance are specified. 

Applicability of solutions. — To illustrate the applicability 
of the generalized solutions, the curves of figure 11 have 
been applied to the previously considered case of the normal 
impact at an initial vertical velocity of 8.86 feet per second 
for comparison with the more exact solution presented in 
figure 4. I n order to make use of the generalized solutions 

3 Although time-history solutions are presented for values of uo' as small as 0.5, it will be 
noted that values of w lma „ <r ma x, v «, and rj tg are not given for values of wo'<1.5. It can be 
seen from the time histories that the characteristics of the solution in the later stages of the im- 
pact change as uo' becomes small; in particular, ui increases and the curve of wi mo , as a func- 
tion of uo' appears to reach a minimum at some value of wo'^1.5. Furthermore, the later 
stages of the solutions greatly stretch out in time and appear to be almost asymptotic in 
character. Several different analytical, numerical, and analogue methods were applied in an 
attempt to study this phase of the problem further but the extremely slow rate of change of 
the variables in this region prevented successful completion of the solutions. 

287846—54 4 


it is first necessary to approximate the tire force-deflection 
characteristics by a simple linear variation. Two such 
linear approximations which might be considered suitable 
for this purpose are shown in figure 12. Linear approxi- 
mation I is a straight line through the origin having a slope 
a=18.5X10 3 pounds per foot (a'=ad= 41.6X10 3 lb). This 
value of a and the other pertinent landing-gear and impact 
parameters result in a value of the initial dimensionless 
velocity parameter u 0 '= 2.57. Linear approximation II is 
a straight line with slope a=21.3X10 3 pounds per foot 
(a'=47.9X 10 3 lb) which does not pass through the origin 
but intersects the displacement axis at a value of 
Z2 (f v _ 0 -j= 0.0508 foot. With this value of a, u 0 '= 2.39. 

Since the solutions of figure 11 have been calculated only 
for integral values of u 0 ' , curves for the foregoing values of 
u 0 ' were graphically interpolated by cross plotting. These 
results were then converted to dimensional values by multi- 
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(d) Relationship between upper-mass acceleration and strut stroke. 
Figure 11. — Continued. 


plying the dimensionless variables by the appropriate 
constants. The results obtained are compared in figure 13 
with the more exact solution presented in figure 4. The 
values based on linear approximation I have been plotted 
exactly as determined from the generalized solutions. The 
results for linear approximation II, however, have been 
displaced relative to the origin of coordinates as indicated 
in the following discussion. 

The assumption of linear approximation II implies that 


the system must move a distance equal to z 2 


( r V° ) 


after 


initial contact (at constant velocity since the wing lift is 
taken equal to the weight) before any finite ground reaction 
can develop. The derivation of the equations of motion, 
on the other hand, assumes that the ground reaction in- 
creases linearly with deflection from the instant of initial 
contact. As a result, the equations of motion do not apply 
until after the system has attained a displacement equal to 


£2 


(?v e ~°y 


which occurs at a time after initial contact 




y 


In other words, the equations of motion 


apply to a coordinate system transformed so that the tire 
force-deflection relationship passes through the origin; that 
is, a coordinate system displaced by z 2(F =0 . relative to 

\ v g ) 


the coordinate system originating at the point of initial 
contact. It therefore follows that the upper-mass and 
lower-mass displacements determined from the generalized 
solutions for the case of linear approximation II must be 
increased by a constant amount equal to z 2(Fv in this 

case 0.0508 foot, and all results must be displaced in 


Z2 (f v =o\ 

time by a constant increment A£= — — -i in this case 


0.0508 

8.86 


=0.0057 second, relative to the instant of initial 


contact. These corrections have been incorporated in plot- 
ting the curves for linear approximation II shown in figure 13. 

As can be seen, the results obtained by application of the 
generalized solutions, particularly by the method employing 
linear approximation II, are in fairly good agreement with the 
more exact solution. The discrepancies which exist are 
attributable to the neglect of the shock-strut preloading and 
springing provided by the air-pressure force, neglect of the 
lower mass, and to differences between the very simple tire 
force-deflection relationships assumed and the exact tire 
characteristics. On the whole, it appears that the general- 
ized results offer a means for rapidly estimating the behavior 
of the landing gear within reasonable limits of accuracy and 
may therefore be useful for preliminary design purposes. 
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(e) Relationship between upper-mass acceleration and upper-mass displacement. 

Figure 11. — Continued. 


SUMMARY OF RESULTS AND CONCLUSIONS 

A theoretical study has been made of the behavior of the 
conventional type of oleo-pneumatic landing gear during the 
process of landing impact. The basic analysis is presented 
in a general form and treats the motions of the landing gear 
prior to and subsequent to the beginning of shock-strut 
deflection. In the first phase of the impact the landing gear 
is treated as a single-degree-of-freedom system in order to 
determine the conditions of motion at the instant of initial 
shock-strut deflection, after which instant the landing gear is 
considered as a system with two degrees of freedom. The 
equations for the two-degree-of-freedom system consider 
such factors as the hydraulic (velocity square) resistance of 
the orifice, the forces due to air compression and internal 
friction in the shock strut, the nonlinear force-deflection 
characteristics of the tire, the wing lift, the inclination of the 
landing gear, and the effects of wheel spin-up drag loads. 

The applicability of the analysis to actual landing gears 
has been investigated for the particular case of a vertical 
landing gear in the absence of drag loads by comparing 
calculated results with experimental drop-test data for corre- 
sponding impact conditions, for both a normal impact and a 
severe impact involving tire bottoming. 

Studies have also been made to determine the effects of 
variations in such parameters as the dynamic force-deflection 


characteristics of the tire, the orifice discharge coefficient, and 
the effective poly tropic exponent for the air-compression 
process, which might not be known accurately in practical 
design problems. 

In addition to the more exact treatment an investigation 
has also been made to determine the extent to which the 
basic equations of motion can be simplified and still yield 
useful results. Generalized solutions of the simplified 
equations obtained are presented for a wide range of landing- 
gear and impact parameters. 

On the basis of the foregoing studies the following con- 
clusions are indicated: 

1. The behavior of the landing gear as calculated from the 
basic equations of motion was found to be in good agreement 
with experimental drop-test data for the case of a vertical 
landing gear in the absence of drag loads, for both a normal 
impact and a severe impact involving tire bottoming. 

2. A study of the effects of variations in the force-deflection 
characteristics of the tire indicates that 

a. In the case of a normal impact without tire bottoming, 
reasonable variations in the force-deflection characteristics 
of the tire have only a relatively small effect on the calculated 
behavior of the landing gear. Approximating the rather 
complicated force-deflection characteristics of the actual tire 
by simplified exponential or linear-segment variations appears 
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o 



Initial dimensionless velocity, 


(f) Variation of maximum upper-mass acceleration and time to reach maximum upper-mass acceleration with initial velocity parameter. 

Figure 11. — Continued. 


to be adequate for practical purposes. Tire hysteresis was 
found to be relatively unimportant. 

b. In the case of a severe impact involving tire bottoming, 
the use of simplified exponential and linear-segment approxi- 
mations to the actual tire force-deflection characteristics 
which neglect the effects of tire bottoming, although adequate 
up to the instant of bottoming, fails to indicate the pro- 
nounced increase in landing-gear load which results from 
bottoming of the tire. The use of exponential or linear- 
segment approximations to the tire characteristics which 
take into account the increased stiffness of the tire that re- 
sults from bottoming, however, yields good results. 

3. A study of the importance of the discharge coefficient 
of the orifice indicates that the magnitude of the discharge 
coefficient has a marked effect on the calculated behavior of 
the landing gear; a decrease in the discharge coefficient (or 
the product of the discharge coefficient and the net orifice 
area) results in an approximately proportional increase in 
the maximum upper-mass acceleration. 

4. A study of the importance of the air-compression process 
in the shock strut indicates that the air springing is of only 
minor significance throughout most of the impact, and that 
variations in the effective polytropic exponent n between the 
isothermal value of 1.0 and the near-adiabatic value of 1.3 
have only a secondary effect on the calculated behavior of 


the landing gear. Even the assumption of constant air 
pressure in the strut equal to the initial pressure (7&=0) 
yields fairly good results, which may be adequate for many 
practical purposes. 

5. An investigation of the extent to which the equations of 
motion for the landing gear can be simplified and still yield 
acceptable calculated results indicates that, for many prac- 
tical purposes, the air-pressure force in the shock strut can 
be completely neglected, the tire force-deflection relationship 
can be assumed to be linear, and the lower or unsprung mass 
can be taken equal to zero. 

6. Generalization of the equations of motion for the 
simplified system described in the preceding paragraph 
shows that the behavior of this system is completely deter- 
mined by the magnitude of one parameter, namely, the 
dimensionless initial-velocity parameter. Solution of these 
generalized equations in terms of dimensionless variables 
permits compact representation of the behavior of the system 
for a wide range of landing-gear and impact parameters, 
which may be useful for rapidly estimating landing-gear 
performance in preliminary design. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 1 , 1952. 
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(g) Variation of maximum upper-mass displacement, maximum lower-mass displacement, and maximum strut stroke with initial velocity parameter. 

Figure 11. — Continued. 
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Figure 11. — Concluded. 
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Figure 12. — Linear approximations to tire force-deflection character- 
istics used in application of generalized solutions. 


(a) Time history of upper-mass acceleration. 

Figure 13. — Comparison of generalized results and more exact solu- 
tion. \ v 0 = 8.8(> feet per second; C<i = 0.9. 
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(b) Time histories of landing-gear displacements and velocities. 
Figure 13. — Continued. 



(c) Time histories of shock-strut stroke and velocity. 
Figure 13. — Concluded. 


APPENDIX A 

NUMERICAL INTEGRATION PROCEDURES 


As previously noted, most of the specific solutions presented 
in this report were obtained with a numerical integration 
procedure, termed the “linear procedure,” which assumes 
changes in the variables to be linear over finite time intervals. 
With this procedure a time interval e=0.001 second was 
used in order to obtain the desired accuracy for the particular 
cases considered. A few of the specific solutions presented 
were obtained by means of a procedure, termed the “quad- 
ratic procedure,” which assumes a quadratic variation of 
displacement with time for successive intervals. This pro- 
cedure, although requiring somewhat more computing time 
per interval, may permit an increase in the interval size for 
a given accuracy, in some cases allowing a reduction in the 
total computing time required. In the case of the more 
exact equations of motion the accuracy of the quadratic pro- 
cedure with a time interval of 0.002 second appears to be 
equal to that of the linear procedure with an interval of 
0.001 second. Although the accuracy naturally decreases 
with increasing interval size, the loss in accuracy for pro- 
portionate increases in interval size appears to be smaller for 
the quadratic than for the linear procedure. In the case of 
the simplified equations of motion reasonably satisfactory 
results were obtained in test computations with the quadratic 
procedure for intervals as large as 0.01 second, whereas the 
linear procedure was considered questionable for intervals 
larger than 0.002 second. 

The generalized solutions presented, because of the 
relatively simple form of the equations of motion, were 
obtained with the well-known Runge-Kutta procedure. 
A study of the allowable interval size resulted in the use of 
an interval A6 = 0.08, which corresponds to a time interval of 
about 0.005 second for the landing gear under consideration. 

UN EAR PROCEDURE 

In this step-by-step procedure the variations in dis- 
placement, velocity, and acceleration are assumed to be 
linear over each finite time interval e. The method, as 
used, involves one stage of iteration. Linear extrapolation 
of the velocity at the end of any interval is used to obtain 
estimated values of velocity and displacement for the next 
interval. These values are then used to calculate values of 
the acceleration in accordance with the equations of motion. 
Integration of the acceleration provides improved values of 
the velocity and, if desired, the displacement and accelera- 
tion. In this procedure all integrations are performed by 
application of the trapezoidal rule. 

The following derivation illustrates the application of the 
linear procedure to the equations of motion for the landing 
gear, which apply subsequent to the beginning of shock-strut 
deflection at time t T . In the example presented internal 
friction forces are neglected in order to simplify the deriva- 
36 


tion. However, the same general procedure can be used if 
these, or other complicating effects, are included in the 
equations. 

For the case under consideration the equations of motion 
(eqs. (16), (17), and (8)) can be written as follows: 

~ 1 £id“2l(ii — i 2 ) 2 +R[l — C(z\ — Z 2 )] ~ n D = 0 (Al) 

4 i ^-A(i 1 -z 2 ) 2 -J5[i-(7(3 1 -2 2 )]-“+rK g (2 2 )-ir 2 =o 

(A 2) 

_ zH — - ‘ Zz-\-Fv (,z$)-\-E=0 (A3) 

g g * 

where 

A 8 pA h 3 

|«| 2(C d A n ycos<p 

B=Pa 0 A a COS<p 


V 0 COSip 

d=k l w-w x 

E=\Y{K l -\) 

Solving equation (A3) for 'i\ gives 

h=[F-Gz 2 -HFy(z,)\\ (A 4) 

g * g 


where 


r=~ (i -K L ) g 


G= 


w 2 

Wi 



Integrating equation (A4) with respect to t between the 
limits l T and t and noting that Zi t =z 2t =z t gives 


Z\— z t -\tFt—G{z 2 


z t ) II J F Vg (z 2 )dr 


(A 5) 


where r= (t — t r ). 

Integrating again and noting that Zi t = z 2r — z T gives 


Zl = (l +G r )(Z T +2rT)+ — Gz 2 — J Fy g (Z 2 )dl 


(A6) 
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Substituting for z x and Z\ in equation (A 2) gives 
f= (a [(1 +G) ( z t -z 2 )+Ft -II j* Fy g (Z 2 ) dr] *+ 

B\l-C^{\+G){Zr+Zrr-Z i n F ^-H^^F Vg { Z2 )dT d^-Fy^+W^ 


(A7) 


The motion of the landing gear subsequent to the 
of shock-strut deflection is determined by means 
by-step solution of equation (A7). This numerical 
yields time histories of the lower-mass motion 
z 2 , z 2} and z 2 , from which the motion variables for 
mass Zi, Zi, and z x can be calculated by means of 
(A4), (A5), and (A6). 

The initial conditions for the step-by-step proc 

Z\ =z 2 „ ='* 


*'n-0= **»-*=*' / 
2 n = 0 = Zr ) 


beginning 
of a step- 
procedure 
variables 
the upper 
equations 

■edure are 


(AS) 




where z T) z T) and z T are the conditions of motion at the 
beginning of shock-strut deflection as determined from the 
solution for the one-degree-of-freedom system. 

Estimated values of the lower-mass velocity at the end of 
the first time increment e following the beginning of shock- 
strut deflection can be obtained from the expression 


H n=l =Zr+eZr 

or, as a first approximation, 

z * nssl =^r 

The corresponding displacement is given by 


(A9) 


Z2„ 


: (Z 2 . + ir) 


(A10) 


After the initial conditions and the conditions at the end 
of the first time increment are established, a step-by-step 
calculation of the motion can be obtained by routine opera- 
tions as indicated by the following general procedure which 
applies at any time T=ne after the beginning of the process. 
The operations indicated are based on integration by appli- 
cation of the trapezoidal rule: 

= + i 2n-2) = ^2»- 1 +| + ^ 2 „_ 2 ) (All) 

« e 2 

^ n =Z2 n _ i +^(Z2„_ , + (Z2„_ } + Z* n _ 2 ) 

(A 12) 

With the estimated values i* and zt the acceleration of 

n n 

the lower mass can be determined by substitution in the 


appropriate integrodifferential equation for the system, 
equation (A7) in the present case. Thus 


Z 2 =.Uz* n , T n ) 


(A 13) 


In equation (A7) the integral expressions can also be evalu- 
ated by application of the trapezoidal rule. For example, 
when F v (z 2 )=mz 2 r , 


rne € 

J o Z 2 r dT~-(z/ + 2z 2l r + . . . 2 Z 2n _ 1 '+Z n r ) 
r(n-1)e 

“Jo + 2 ,0 


(A 14) 


and 


• (n-l)e 


f*ne (*ne f*(n— l)e ( 

Z 2 d t dr ~ I I z 2 dr dr T 

Jo Jo Jo Jo 

f / r(n-l)e r ne \ 

id ’> dr +i ! < d ') 


(A 15) 


An improved value for the velocity is obtained from the 
expression 


22 B =2 2„_,+f + 


(A 16) 


This value is used in the calculation of the estimated velocity 


~*n + l 


and displacement z* n+1 for the next interval. 


If desired, improved values of the displacement and 
acceleration for the nth interval subsequent to the beginning 
of shock-strut deflection can be obtained as follows: 

Z2„=Z 2ri _,+^ (Z2 n _ 1 +Z 2 „) 



€ 2 

= 22„_, + « 2 2„_,+ J ( Z 2 b _ 1 + 2 2 „) 

(A 17) 

and 




**. = R*V 2 V T ») 

(A 18) 


where j(z 2n , z 2n , r„) is an appropriate equation for the system, 
such as equation (A7). 

With the values of z 2 , z 2 , and z 2 , the motion variables for 

-tv -tv -tv 

the upper mass z ln , z lfi , and z lfj can be calculated separately 
from equations (A4), (A5), and (A6), as previously noted. 

In setting up the numerical procedure used in obtaining the 
solutions presented in this report, an evaluation of the errors 
introduced by the procedure indicated that it would not be 
necessary to calculate the improved values of the displace- 
ment z 2n (eq. (A17)) or the acceleration z 2n (eq. (A18)). How- 
ever, improved values of the velocity z-? n were calculated by 
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means of equation (Al6) for the purpose of determining 
estimated values of the velocity z 2 and the displacement z 2 
(eqs. (All) and (Al2)) for the increment immediately 
following. 

In order to illustrate the application of the method, a 
tabular computing procedure for the solution of the system 
represented by equations (Al), (A2), and (A3) is presented in 
table I. 

QUADRATIC PROCEDURE 

In this step-by-step procedure a quadratic variation of 
displacement is assumed over successive equal finite time 
intervals for the purpose of extrapolating values of the 


motion variables from one interval to the next. With this 
assumption the displacement variation over two successive 
equal time intervals is completely determined by the three 
values of displacement at the beginning and end of each of 
the two intervals. By writing the quadratic variation in 
difference form, the velocity and acceleration at the midpoint 
of the double interval can be expressed in terms of the three 
displacement values previously mentioned. Substituting for 
the velocity and acceleration in the differential equations for 
the system yields difference equations of motion in terms of 
successive displacement values which can be evaluated 
interval by interval. 


TABLE I 

LINEAR PROCEDURE 


Row 

Quantity 

Equation 

Procedure f 

(0) 








© 

Z*2 n 

*2.— l + f (^._2+22„-l ) 

®„ 

© 



©p + C® p 

© 

F ^n) 


Determined from tire force- 
deflection characteristics. 

© 

J >'«* 

Equation (A 14) 

®»+| [®+©p] 

© 

rw-yo'* 

Equation (A 15) 

©p+|[©+®p] 

© 

K 

g 

Equation (A7) 

Given by equation (A7). 

© 

K 

^„.,+ 2 (&„_, + %„) 

©p+2 [©+©«>]£ 

© 

2 (*!„+ 

Z 2n +|(Z2 n _, + 22 n ) 

©+{ [©+©pk 

© 

Z - n-H 

^>, + 1 (^’„-i + 22„) 

©+| [®+©pk 

@ 

.. 

fu 

g 

Equation (A4) 

Given by equation (A4). 

® 


Equation (A5) 

Given by equation (A5). 

© 

Zi n 

Equation (AO) 

Given by equation (AO). 


t O p denotes value for previous time interval. 
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The following derivation shows how the procedure can be 
applied to the determination of the behavior of the landing 
gear subsequent to the beginning of shock-strut deflection at 
time t T . In order to simplify the derivation, internal friction 
forces are again neglected in setting-up the equations of 
motion. 

The assumption of a quadratic variation of displacement 
with time (constant acceleration) over two successive inter- 
vals, each of duration e, permits expressing the velocity and 
acceleration at the midpoints of the double interval (see 
sketch) in terms of the displacement values at the beginning, 
midpoint, and end of the double interval by the equations 
(see ref. 5, p. 16): 

in= Zn - +l ~ e Zn ~ l (A 19) 

and 


Z n ~ 


z n+1 — 2z n + z n _! 


(A20) 



where z n , z n , and z n are the velocity, acceleration, and displacement at the end of the nth interval ( r=ne ) after the 
beginning of shock-strut deformation and z n _ x and z n +i are the displacements at the end of intervals n— 1 and n + 1, 
respectively. 

Substituting the difference relations for z i, i 2 , Zi, and z 2 into equations (Al) and (A3) permits writing the equations 
of motion for the landing gear in difference form as follows: 


Wi, 

ge 2{Zi n + 1 


2Z l* + Zl »-l)+ 472 ( Z l* +I -Zl„-l — S 2* + 1 + 2 2 K _,) 2 + £[ 1 -C(Zl n — Z2„)] “+£> = 0 


and 

z i»+i =2z >„- z i„— i- Gr ( z 2»+i-2 Z 2 a +z 2n _ 1 )-He 2 [F Vt (z t „')+E] 

where the constants are as defined in the previous section. 

Substituting for z ln+1 in equation (A21) gives 

z ’n+ i=% i+ ^V 2 [2WVTr 2 _ A 4 WSW&A WK^WTW^A TT*(4 W r 1 *«, + i+^T. + i)] 


(A21) 

(A22) 


(A23) 


where 

a n+ i=2W 2 Z2 n — W 2 Z2 n _ 1 —ge 2 [F v g (z 2 ji )+E] 

(3 n + 1 = 2 IT 2 Z 2 n +(W 1 -—W 2 )z 2n _ l -\- 2 Wi(zi n —Zi n '_ 1 )—ge 2 [F Vg (z 2n ')-t-E] 

and 

Equations (A.22) and (A23) are essentially extrapolation formulas which permit the determination of values for the 
upper-mass and lower-mass displacements to come from the values of displacement already calculated. These equations 
thus permit step-by-step calculation of t lie displacements as the impact progresses, starting with t lie initial conditions, 
from which the upper-mass and lower-mass velocities and accelerations can be determined by means of equations (A 19) 
and (A20). 

Since the calculation of the displacements Z\ and z 2 at any instant by means of equations (A22) and (A23) requires 
values for the displacements at two previous instants, the routine application of these equations can begin only at the 
end of the second interval (r=2e) following the beginning of shock-strut deflection. Before the displacements at the end 
of the second interval can be calculated, however, it is necessary to determine the displacements at the end of the first 
interval. These values can be obtained from the conditions of motion at the instant of initial shock-strut deflection by 
applying equations (Al9) and (A20) to the instant t=t T . 
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At the instant of initial shock-strut deflection 


Zl n = 0 ^ 2 n = 0 ^ T 

Z'n = 0 =i *n = 0 =i W 
2l n=0 ~ 2 n = 0 ~V 


(A 24) 


Application of the difference equations (A19) and (A20) 
to the instant t=t T (that is, ?i=0) gives the following 
equations: 

Zn = 1 ^ n = — 1 i 


2t = - 


Z T = ~ 


Zn-\~ 2Zr-\-Z n= - 


S 


(A 2 5) 


Since the landing gear is considered as a one-degree-of- 
freedom system from initial contact up to the instant t=t T , 


the foregoing application of the difference equations results 
in identical values for the upper-mass displacement and 
lower-mass displacement at t lie end of the first interval. 
Simultaneous solution of equations (A25) gives the following 
expression for the displacement at the end of the first interval: 

22 B= 1 =2. n = 1 = Zr+eir+|- 2 r (A26) 

With the values for e T and z n= i, equations (A22) and 
(A23) permit the step-by-step calculation of the upper-mass 
and lower-mass displacements subsequent to the first interval 
following the beginning of shock-strut deflection. The 
corresponding velocities and accelerations of the upper and 
lower masses can be determined from the calculated displace- 
ments by means of equations (A19) and (A20), as previously 
noted. 

A tabular computing procedure illustrating the application 
of the method is presented in table II. 


TABLE II 

QUADRATIC PROCEDURE 


Row 

Quantity 

Equation 

Procedure t 


T 





® 



©p 

© 

2| n 


®p 

© 

*„+, 

Equation (A23) 

Given by equation (A23). 

© 

2l n+l 

Equation (A22) 

Given by equation (A22). 

© 

Zo 

n 

Z *n+ \~ Z ' 2 n-\ 
2 6 

®-@p 

2 e 

© 

Z2 n 

Z 2„+l“ 22 ->n+^n-l 

€ 2 

®-2®+ ® p 

e 2 

© 


2 '„ + .- 2 ' n -i 

2 e 

®-®p 

26 

© 

K 

2 l» +I - 22 '»+ 2 *»-l 

€ 2 

®— 2© + ©„ 

e 2 

1 


f O p denotes value for previous time interval. 
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RUNGE-KUTTA PROCEDURE 

In this step-by-step procedure the differences in the de- 
pendent variables over any given interval of the independent 
variable are calculated from a definite set of formulas, the 
same set of formulas being used for all increments. Thus the 
values of the variables at the end of any given interval are 
completely determined by the values at the end of the pre- 
ceding interval. Unfortunately, however, unless the equa- 
tions to be integrated are relatively simple, the method can 
become quite lengthy. 

The following derivation illustrates the application of the 
Runge-Kutta method to the generalized equations of motion 
(eqs. (21)) for the simplified system considered in the section 
on generalized results. Since these equations can be readily 
reduced to the first order, they can be integrated by the 
step-by-step application of the general equations given on 
pages 301 and 302 of reference 6 for first-order simultaneous 
differential equations. 

The generalized equations for the simplified system pre- 
viously discussed (eqs. (21)) are 

(ui— u 2 ') 2 +Ui"=0 

(ui — u 2 ') 2 —u 2 =0 


U\ f ~\- / u 2 -— 0 


Inasmuch as any two of these equations are sufficient to 
describe the behavior of the system, only the last two equa- 
tions are employed in this procedure. These equations can 
be reduced to a first-order system by introducing the new 
variable 

w=Ui (A27) 

so that 

w' —Ui" (A 2 8) 

and the equations of motion become 


(w—u 2 / ) 2 —u 2 =0) 

\ (A 2 9) 

w'+u 2 =0 ) 

Solving equations (A29) for u 2 ' and w', respectively, gives 

u/=w-^ (A30) 

v/ = —u 2 (A31) 


and (A31) gives 

Aui = Ui n — Ui nl =— (ki~\~2k 2 + 2£ 3 — |- k 4 ) 
A w=w n — w w _ 1 = g (/i + 2/2U 2/3 +/ 4 ) 




► (A 3 2) 


where 


Au 2 =u 2n — w 2n _ 1 =- (m 1 +2m 2 -f 2m 3 +m 4 ) 


&i=w„_iA0 


fc 2 =(w»-i+|) A0 

k 3=(w„_i+!) A0 

A-4=(w n _,-M 3 )A0 


h=—u 2n _ l Ad 

^2=— («2„_,+^ i ) A 0 

l3=—(u 2 „_ !+-y) A0 
h= — (u 2n _ l + m 3 ) A0 

m l = (w„_ l —^/u 2n _ l ) A0=m 2 ' _ 1 A0 

m 3 =j^(w n _ i +|)— -^/u 2;( _ 1 +^J A0 

m 4 =[(w n -i+W— -V w 2 n -i+ m 3 ] Ae 

With this procedure, Ui, w, and u 2 can be calculated in 
step-by-step fashion from the values for the preceding inter- 
val, the procedure beginning with the initial conditions. 
From these values, Ui, U\\ and u 2 can be calculated by 
means of equations (A27), (A28), and (A30), respectively. 


Applying the general procedure presented in the reference 
previously cited to the simultaneous equations (A27), (A30), 


APPENDIX B 


SOURCE OF EXPERIMENTAL DATA 


Following is a brief description of the apparatus and test 
specimen used in obtaining the experimental data presented 
in this report. 

EQUIPMENT 

The basic piece of equipment employed in the tests is the 
carriage of the Langley impact basin (ref. 7) which provides 
means for effecting the controlled descent of the test speci- 
men. In these tests the impact-basin carriage was used in 
much the same manner as a conventional stationary landing- 
gear test jig (see ref. 8). In order to simulate mechanically 
the wing lift forces which sustain an airplane during landing 
the pneumatic cylinder and cam system incorporated in the 
carriage was used to apply a constant lift force to the dropping 
mass and landing gear during impact. The lift force in these 
tests was equal to the total dropping weight of 2,542 pounds. 


TEST SPECIMEN 


The landing gear used in the tests was originally designed 
for a small military training airplane having a gross weight of 
approximately 5,000 pounds. The gear is of conventional 
cantilever construction and incorporates a standard type of 
oleo-pneumatic shock strut. The wheel is fitted with a 27- 
inch type I (smooth-contour) tire, inflated to 32 pounds per 
square inch. The weight of the landing gear is 150 pounds. 
The weight of the lower mass (unsprung weight) is 131 
pounds. 

In the present investigation the gear was somewhat modi- 
fied in that the metering pin was removed and the original 
orifice plate was replaced with one having a smaller orifice 
diameter. Figure 14 shows the internal arrangement of the 
shock strut and presents details of the orifice. Other perti- 
nent dimensions are presented in table III. The strut was 
filled with specification AN-W-0-366B hydraulic fluid. 
The inflation pressure with the strut fully extended was 43.5 
pounds per square inch. In these tests the landing gear was 
mounted with the shock-strut axis vertical. Figure 15 
is a photograph of the landing gear installed for testing. 


TABLE III 

IMPORTANT CHARACTERISTICS OF LANDING GEAR 
USED IN TESTS 


A a , sq ft 
A h , sq ft 
A 0 , sq ft 

Vo, CU ft _ __ 

Pa 0 , lb/sq ft 

h, ft-. 

/, ft 
W u lb 

W 2 , lb 


0. 05761 
0. 04708 
0. 0005585 
0. 03545 
6, 264 
0. 5521 
2. 22604 
2, 411 
131 


INSTRUMENTATION 

A variety of time-history instrumentation was used during 
the tests. The vertieal acceleration of the upper mass was 
measured by means of an oil-damped electrical strain-gage 
accelerometer having a range of ±85 and a natural frequency 
of 85 cycles per second. A low-frequency (16.5 cycles per 
second) NACA air-damped optical-recording accelerometer, 
having a range of — 1 g to 6 g, was used as a stand-by instru- 
ment and as a check against the strain-gage accelerometer. 
Another oil-damped strain-gage accelerometer, having a 
range of ±12® and a natural frequency of 260 cycles per 
second, was used to determine the vertical acceleration of the 
lower mass. The vertical displacement of the lower mass 
(tire deflection) and the shock-strut stroke were measured 
separately by means of variable-resistance slide-wire poten- 
tiometers. The vertical displacement of the upper mass was 
determined by addition of the strut-stroke and tire-deflection 
measurements. The vertical velocity of the landing gear at 
the instant of ground contact was determined from the output 
of an elemental electromagnetic voltage generator. A time 
history of the vertical velocity of the upper mass was ob- 
tained by mechanically integrating the vertical acceleration 
of the upper mass subsequent to the instant of ground con- 
tact. Electrical differentiation of the current output of the 
strut-stroke circuit provided time-history measurements 
of the shock-strut telescoping velocity. The instant of 
ground contact was determined by means of a micro- 
switch, recessed into the ground platform, which closed 
a circuit as long as the tire was in contact with the platform. 

The electrical output of the instruments was recorded on a 
14-channel oscillograph. The galvanometers were damped 
to approximately 0.7 critical damping and had natural fre- 
quencies high enough to produce virtually uniform response 
up to frequencies commensurate with those of the measuring 
instrumentation. A typical oscillograph record is shown in 
figure 16. 

It is believed that the measurements obtained in the tests 


are accurate within the following limits: 

Measurement Accuracy 

Upper-mass acceleration, g ±0. 2 

Force on upper mass, lb ± 500 

Lower-mass acceleration, g ±0- 3 

Vertical velocity at ground contact, fps _ ±0. 1 

Upper-mass velocity during impact, fps_ ±0. 5 

Upper-mass displacement, ft ±0.05 

Lower-mass displacement, ft ±0. 03 

Shock-strut stroke, ft ±0.03 

Shock-strut telescoping velocity, fps ±0. 5 

Time after contact, sec ±0. 003 
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Filler plug 



Figure 14. — Shock strut of landing gear tested at Langley impact basin. 
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Figure 15. — View of landing gear and instrumentation. 
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Figure 16. — Typical oscillograph record obtained during test in Langley impact basin. 
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